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A possible solution to the observed baryon asymmetry in the universe is described, based on the 
physics of the standard model of electroweak interactions. At temperatures high enough electroweak 
physics provides violation of baryon number, while C and CP symmetries are not exactly conserved, 
although in the context of the minimal electroweak model with one Higgs doublet the rate of CP 
violation is not sufficient enough to generate the observed asymmetry. The condition that the 
universe must be out of thermal equilibrium requires the electroweak phase transition (EWPT) to 
be first order. The dynamics of the phase transition in the minimal model is investigated through 
the effective potential, which is calculated at the one loop order. Finite temperature effects on the 
effective potential are treated numerically and within the high temperature approximation, which is 
found to be in good agreement with the exact calculation. At the one loop level the phase transition 
was found to be of the first order, while the strength of the transition depends on the unknown 
parameters of the theory which are the Higgs boson and top quark masses [l| . 
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"Why does the whole world have < <p >= +v? Why 
doesn't it have < <fi >= —v somewhere? Suppose that God 
created the universe in the state < <j) >— and then the 
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universe discovered that it could lower its energy; Where 
it puts its energy is none of my business, but it gets rid 
of it - gives it back to God or something; " 

R. Feynman 3] 

I. THE BARYOGENESIS PROBLEM 
A. Overview 

Recently the connection between particle physics and 
cosmology has received much attention. The classical 
cosmological model of the expanding universe provides 
a powerful framework where particle physics can test its 
predictions about matter genesis. The expansion of the 
universe can be considered as an enormous particle ac- 
celerator, which although ran billions of years ago in the 
far past, can still be used to check the validity of theo- 
ries concerning the elementary particles. Cosmology on 
the other hand can use the predictions of particle physics 
about the nature and behaviour of the elementary par- 
ticles in order to cure unsolved up to now cosmological 
problems which are involved in the theories concerned 
with the evolution of the universe. 

The baryogenesis problem, that is the observed excess 
of matter over antimatter today in the universe, is one 
of these problems. As pointed out by Sakharov in 1967 
0, the universe began in a baryon symmetric state but 
particle interactions produced a net asymmetry. He pos- 
tulated three conditions to be satisfied in order to explain 
the observed baryon excess. These are: (a) baryon num- 
ber non-conservation, (b) C and CP violation and (c) 
that the universe must be out of thermal equilibrium. 
Our aim is to explore how these conditions are satisfied 
in the framework of the standard model of electroweak 
interactions. We will focus our attention on the third 
condition and provide a careful study of the electroweak 
phase transition. Our main goal is to examine the order 
of this transition when temperature is introduced in the 
model via the so-called one loop approximation. 

This work is organized as follows: In order to estab- 
lish the basic principles for the baryogenesis problem the 
remainder of this section is devoted to a brief account 
on the development of the big bang model, the standard 
cosmological model, from its beginning up to the present 
days. The most important stages in the development of 
the big bang model are given in Section lTrH In Scction lTCl 
we present how the baryogenesis problem appears in the 
framework of the standard big bang cosmology. The ne- 
cessity to satisfy the Sakharov conditions and relevant 
comments are outlined in Section ll PI A recent develop- 
ment in the standard hot big bang model, the inflation 
model, is discussed briefly in Section ITeI 

According to recent investigations it is possible that 
the necessary conditions for the solution of the baryo- 
genesis problem can be satisfied in the standard model 
of the electroweak interactions, so we find very useful to 
provide an introduction to the standard model of these 



interactions and its connection with cosmology. This is 
the subject of Sectional where phase transitions in gauge 
theories are also discussed briefly through well known ex- 
amples. 

The standard electroweak model satisfies the first two 
conditions, so we sketch in Section IlIII how baryon num- 
ber and CP symmetry are violated in the model. In the 
standard model, due to the non-trivial structure of the 
SU(2) gauge vacuum and as a consequence of anomalous 
processes, baryon number is not conserved. On the other 
hand CP symmetry is violated in the standard model 
because of relative phases between the electroweak gauge 
interactions and the Higgs interactions of the quarks. 

The third condition is satisfied if the electroweak phase 
transition (EWPT) is of the first order. The basic tool 
for the investigation of the electroweak phase transition 
is the finite temperature effective potential, the quantity 
which has the meaning of the free energy density of the 
system under consideration, and the relevant theory is 
given in some detail in Section ll VI 

The evolution and the order of the electroweak phase 
transition will be determined in the so-called one loop 
approximation in Section [V] Finite temperature effects 
are treated exactly and within the high temperature lim- 
iting case. 

The last section is devoted to a discussion concerning 
our investigation of the electroweak phase transition in 
connection with the baryogenesis problem and we also 
present there a summary of our results and conclusions. 



B. The Hot Big Bang 

Modern cosmological theories started to develop within 
the framework of Einstein's general theory of relativity 
early this century. The development of particle physics 
in the recent years has infused new ideas in cosmology, 
which can be applied to the study of the earliest mo- 
ments of the universe and the evolution to its present 
state. The most accepted up to date cosmological theory 
is the hot big bang model. This theory is based upon 
general relativity and the Friedmann model for an ex- 
panding universe. According to this model the universe 
starts its life from a state of enormous matter density 
and temperature and after its expansion results in what 
we observe today. The evolution and the most important 
stages in the development of this cosmological theory up 
to the present days have been explored in an elegant way 
by many authors [HIS0- 

A detailed discussion on the formulation of this theory 
is given by Weinberg Q and Kolb and Turner , while a 
brief account can be found in Lindc The beginning 
of the theory and the first important stage goes back in 
1922 when Friedmann based on Einstein's cosmological 
principle, which is the hypothesis that the universe is spa- 
tially homogeneous and isotropic, and created a model for 
an expanding universe. A few years later, Hubble's and 
Slipher's observations of the galactic redshift supported 
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experimentally the idea of the expansion from an initial 
state and very soon this theory became the cornerstone 
of modern cosmology. 

In Friedmann's model the radius of the universe R = 
R(t), or more explicitly its scale factor, is a function of 
time and has an evolution which is described by the equa- 
tions of general relativity. Einstein's equations imply an 
evolution for the scale factor given by the first order dif- 
ferential equation 

R + k= ^ pR 2. (1) 

The constant k in the above equation describes the cur- 
vature of space and takes the values k = — 1, 0, 1, which 
correspond to open, flat and close universe models respec- 
tively and G is the gravitational constant. In addition we 
have an energy conservation law which can be expressed 
by 

pR 3 + 3(p + P )R 2 R = , (2) 

where p is the energy density of matter in the universe 
and p its pressure. In order to solve these equations and 
find out how the universe evolves in time, we need a state 
equation which relates the energy density of the universe 
to its pressure. We may assume p — top so the above 
equations give the result 

p ~ iT 3(1+w) . (3) 

We can consider two possibilities, in order to interpret 
this solution. If we suppose that the universe was filled 
with nonrclativistic cold matter with negligible pressure 
p = 0, we find that p ~ R~ 3 . If the universe was a 
hot ultrarelativistic gas of noninteracting particles with 
p = p/3, we then find that p ~ R~ A . When R is small 
one finds that for nonrelativistic cold matter the radius 
of the universe evolves in time as R ~ i 2 / 3 , while in the 
case of the hot ultrarelativistic gas it goes as R ~ t 1 ' 2 . 
Thus regardless the model used, as the time t tends to 
zero, the scale factor vanishes and the density of matter 
becomes infinite. The point t = is known as the initial 
cosmological singularity. It is the initial point when the 
universe starts its life and then expands to what it is 
today. The expansion rate of the universe is given by 
the quantity H = Rj R which is known as the Hubble 
constant, although it is not really a constant but it varies 
with time. 

Up to the mid-sixties it was not clear if the universe 
had started its life from a hot or a cold initial state. The 
new era in cosmology opened when Penzias and Wilson 
discovered the 2. 7° if microwave background radiation in 
1964-65, which had been predicted by the hot universe 
theory. If we suppose that the universe expanded adi- 
abatically, then during the expansion the quantity RT 
remained approximately constant, and the temperature 
of the universe dropped off as T ~ 1 (£) . The radi- 
ation which was detected by Penzias and Wilson were 
the relic photons of the initially hot photon gas which 



cooled down during the expansion of the universe. This 
was the second stage in the development of modern cos- 
mology and was decisive in the establishment of the hot 
big bang model as the standard cosmological model. 



C. Baryon Asymmetry in the Universe 

Although the standard hot big bang model is suc- 
cessful in accounting for Hubble expansion, the residual 
microwave background radiation and the abundances of 
light atomic nuclei there remain some serious difficulties 
and one of these open questions is the observed baryon 
asymmetry in the universe today. The essence of the 
baryon asymmetry problem is to understand why in the 
observable part of the universe the density of baryons is 
many orders greater than the density of antibaryons and 
why on the other hand the density of baryons is much less 
than the density of photons The baryon asymmetry 
can be quantified by a quantity known as the baryon 
number B of the universe and is defined as the ratio of 
the baryon number density ub to the entropy density s. 
The baryon number density ns is defined as the number 
density of baryons rib minus the number density of an- 
tibaryons 7i£, so ns — rib — Ttr. The baryon number of 
the universe today, is estimated in the range |9j, [llj [See 
however in Section IV 1 (Jj 

B = — ~ (6 - 10) x ltr 11 . (4) 

s 

There is a strong evidence of the asymmetry since it is 
well known that the earth as well as all the objects on it 
consist of matter rather than antimatter. It is an experi- 
mental fact that there is no antimatter on the earth and 
the outer space. Measurements of the cosmic rays emit- 
ted from the sun have proven the absence of antimatter 
in the solar system. Only an insignificant positron flux 
is contained in solar cosmic rays and it is mainly due to 
interactions of the primary solar cosmic rays with matter. 

Astronomical observations show that at least our 
galaxy consists of usual matter, while an inconsiderable 
amount of antimatter is observed and does not seem to 
indicate the presence of antimatter in the galaxy. For 
clusters of galaxies and their intergalactic gas one can 
derive limits for the amount of antimatter present from 
the 7 flow expected from the 7r° decay, since matter- 
antimatter annihilation would produce 7r° mesons and 
those decay into 7-rays. These hard photons have never 
been observed. It is not definite if there are in our uni- 
verse islands of antimatter separated by empty space of 
regions containing normal matter. Electromagnetic ra- 
diation which is the main source of informations from 
the outer space cannot give a signal about its matter or 
antimatter origin. 
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D. Three Cornerstones of Baryogenesis 

In Friedman's cosmology the excess of matter over an- 
timatter is considered as one of the initial conditions 
and the baryonic asymmetry as one of the fundamen- 
tal cosmological constants. This philosophy changed in 
early seventies when Sakharov and others suggested that 
the observed asymmetry could be generated dynamically 
even though the universe had started in an initially sym- 
metric state. Sakharov postulated three necessary con- 
ditions for this to happen: 

(a) There must be baryon number violating 
processes. 

(b) Breaking of charge (C) and combined 
(CP) symmetries. 

(c) These processes must go out of equilib- 
rium sometime during the history of the 
universe. 

The first condition is obvious since it is clear that the 
baryon number must be violated if the universe started 
its life in a baryon symmetric state and then evolved into 
an asymmetric one. Many models predict baryonic asym- 
metry |12| . although baryon number violation has never 
been verified experimentally. In Grand Unified Theories 
(GUT) baryon number violation occurs in the equations 
of motion of the theory. It proceeds with the exchange of 
very massive particles and takes place at an energy scale 
of the order of magnitude about 10 16 GeV. Supersym- 
metric models also predict baryon non conserving pro- 
cesses in the energy range 10 16 GeV- 10 2 GeV depending 
upon the concrete model. 

In order to understand the necessity of the second con- 
dition we note that the baryon number of a state is odd 
under C and combined CP symmetries. That is, the 
baryon number of a state changes sign under charge con- 
jugation C and charge conjugation combined with parity 
CP. Thus if a state is either C or CP invariant, then 
the baryon number of this state should be zero. If the 
universe begins its life with equal amounts of matter and 
antimatter, and without a preferred direction, then its 
initial state is both C and CP invariant. Unless both C 
and CP are violated, the universe will remain C and CP 
invariant as it evolves to its final state, and thus cannot 
develop a net baryon number even if the baryon number 
is not conserved. Therefore both C and CP violations are 
needed if a net excess of matter over antimatter is to be 
produced. In contrast to baryon number violation, which 
has only been predicted in various theoretical models but 
with no experimental evidence for it, C and CP viola- 
tion have been verified experimentally in neutral Kaon 
decays. 

The necessity to satisfy the last Sakharov condition, 
that the universe must be out of thermal equilibrium, is 
very important, since if an initially baryonic symmetric 
universe is in thermal equilibrium, then particle number 



densities are given by 

1 

neq ~ eM(E-»)/T]±l ' (5) 

where E is the particle energy and /i its chemical po- 
tential. If the charge is not conserved the corresponding 
chemical potential vanishes in equilibrium. CPT invari- 
ance guarantees that particles and antiparticles have the 
same mass, therefore their number remains equal during 
the expansion and so no asymmetry occurs, regardless of 
B, C, or CP violating processes. 



E. The Inflation Model 

The development of gauge theories and their influence 
to the cosmological questions led to the development of 
a new version of the standard big bang model. The in- 
flation model which was introduced by Guth in 1981 0] 
came to cure some of the pathological problems of the hot 
universe theory. According to this model, soon after the 
beginning of the big bang the expansion rate of the uni- 
verse increases exponentially, and due to this expansion 
the universe supercools and all the physical processes are 
interrupted. Then the universe reheats again and contin- 
ues its expansion with a slower rate. 

A detailed analysis of the inflation model is given by 
Linde |ToL IT^j . while important papers relevant to this 
model (and not only) are collected by Abbott and Pi [lflj . 
Although the inflation model is not the subject of this 
work, we should note here that it has already solved many 
of the problems that the hot big bang model possesses. 
On the other hand the inflation model imposes the con- 
straint that the baryon generation should have occurred 
after the exponential expansion, since the creation of a 
large amount of entropy during the inflation epoch would 
dilute any pre-existing baryon asymmetry. If one be- 
lieves the inflation model then baryogenesis should hap- 
pen at rather lower scale energies than those of the GUT's 
scale. This is especially gratifying for present day inves- 
tigations, since our technology makes only electroweak 
energy scales, not GUT's scales, accessible to experimen- 
tal tests. 



II. GAUGE THEORIES AND COSMOLOGY 

A. Introduction 

The experimental success of the model of unified elec- 
troweak interactions and the development of unified 
gauge theories of all the interactions at the beginning 
of the seventies opened a new era for the theories of the 
elementary particles. The creation and evolution of these 
theories have their roots on symmetry principles. Gen- 
eralization of the ideas about symmetries generated the 
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gauge principle and the present belief that all particle in- 
teractions are dictated by the so-called local gauge sym- 
metries. This principle arises from the requirement that 
all quantities which are conserved, are conserved locally 
and not merely globally. 

Gauge invariant Lagrangians describing particle inter- 
actions remain invariant under a set of local gauge trans- 
formations. But invariance of the Lagrangian demands 
for the particles to be massless. The particles become 
massive due to the appearance of constant classical scalar 
fields <p over all space, through the so-called mechanism 
of symmetry breaking. 

The rapid development of particle physics not only led 
to a better understanding of particle interactions, but 
also to a significant progress in the theories of superdense 
matter, matter more than 80 orders of magnitude denser 
than the nuclear matter. It was Kirzhnits and Linde 0] 
and very soon others 0, who showed that the scalar 
fields which are responsible for the breaking of symme- 
try, could disappear at high enough temperatures. This 
means that at temperatures high enough phase transi- 
tions take place and after they have been completed, the 
symmetry is restored. 

These theories are applicable in cosmology since su- 
perdense matter at very high energies and temperatures 
should appear at very early stages of the evolution of the 
universe. According to the standard hot big bang model 
for an expanding universe, as the universe cools down 
these phase transitions could take place. Cosmological 
problems, as the baryon excess for example, could find 
an answer in the framework of gauge theories. Recent 
development in this area suggest that the baryogenesis 
problem could find its answer in the standard model of 
electroweak interactions. This model possesses baryon 
number violation as it was proved by t'Hooft [T^. |20| . 
although the violation rate at zero temperatures seems 
to be too small. Violation of the CP symmetry is pro- 
vided by the Cabbibo-Kobayashi-Maskawa mixing be- 
tween the quark families. The condition that the sys- 
tem must be out of equilibrium demands the electroweak 
phase transition to be of the first order. 

The remainder of this section is devoted on the formu- 
lation of the electroweak theory and an introduction to 
the theory of phase transitions in gauge theories. In the 
next section we discuss the formulation of the standard 
electroweak model as a gauge theory. The mechanism of 
symmetry breaking and the generation of particle masses 
are presented in Sections 111 CI and 111 D\ respectively. The 
idea of symmetry restoration at high temperatures and a 
brief discussion of the cosmological phase transitions are 
contained in the last two Sections III El and III Fl . 



B. Electroweak Gauge Theory 

The successful formulation of Quantum Electrodynam- 
ics (QED) as a gauge theory led to the idea of extending 
the gauge principle to the description of other interac- 



tions. Constructing a gauge theory we demand the La- 
grangian to be locally invariant under a group of inter- 
nal symmetries. Thus we introduce into the Lagrangian 
a number of vector fields equal to the number of the 
symmetry group generators. The gauge bosons self- 
interactions as well as the way they couple to matter 
fields is completely determined by the gauge couplings. 
The QED Lagrangian for example is invariant under a set 
of local U(l) transformations with generator Q, the elec- 
tric charge operator. Invariant interactions are mediated 
by one vector field, the photon field. In Quantum Chro- 
modynamics (QCD), the gauge theory of strong interac- 
tions, we introduce eight vector fields (the gluon fields) 
into the Lagrangian and the theory is based on the group 
SU(3), which has eight generators the Gell-Mann matri- 
ces. 

The Glashow-Salam-Weinberg model which unifies 
the weak and electromagnetic interactions, known as the 
standard model of electroweak interactions, is based on 
the non-Abelian group SU(2)l x U(1)y and has four 
generators. The weak interactions are mediated by three 
vector bosons A" a = 1,2,3 which are accommodated 
in the adjoint representation of an SU(2) group denoted 
by SU(2)l with generators T a — r /2, where the r OJ 
a = 1, 2, 3 are the Pauli matrices and the T a are referred 
as the weak isospin generators. The subscript L signifies 
that the A° couple in a parity violating way to the left 
handed parts of the lepton matter fields, which transform 
as a doublet again under the same SU(2)l. 

To accomodate the electromagnetic interactions we 
need another gauge boson and another group de- 
noted by U(1)y- The generator of the U(l)y group is 
Y"/2, where Y is the weak hypercharge defined through 
the so-called Gcll-Mann-Nishijima relation as 

Q = T 3 + Y/2 (6) 

and T3 is the third component of the weak isospin. The 
gauge coupling for the weak isospin group SU(2)l is de- 
noted by g and that for weak hypercharge group U(l)y 
as g'/2. The kinetic terms into the Lagrangian for the 
two gauge fields are 

C — --F a W"_1r rm" (7) 

'-'gauge— ^ ^ J - ) [j t v- LJ j \> ) 

where the field strength tensor for the SU (2) L gauge field 
is given by 

F; v = d»A a v - d v Al + ge^A^Al , (8) 

while the field strength tensor of the U(l)y gauge field 
has the form 

= d„B v - d v B„ . (9) 

The physical particles which mediate the weak and elec- 
tromagnetic interactions, the two charged W + and W~, 
the neutral Z° and the photon 7 result as linear combi- 
nations of the fields A° and J3 M . 
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The fermionic matter fields are present in the model 
grouped in three families or generations which include the 
e, v e , u, d known as first generation and other two repli- 
cations the fjt, v^, c, s family and the r, v T , t, b family. 
Violation of the parity is introduced into the model by 
putting the left-handed and right-handed fermions into 
different group representations. All the left-handed fields 
are taken to transform as doublets under SU(2)l, while 
the right-handed fermions transform as singlets. 

The left-handed quarks are introduced into the model 
grouped into SU (2) doublets of the form 

Qu = ( u d ,) Qc=( c s ,) Q*=(v) • ("J) 

The prime means that these are weak eigenstates, which 
are constructed as linear combinations of the mass eigen- 
states. The matrix which connects the two sets of eigen- 
states is referred as the Cabbibo-Kobayashi-Maskawa 
(CKM) matrix. By convention the charge 2/3 quarks are 
unmixed, but the eigenstates of the charge —1/3 quarks 
are related through the CKM matrix as 




This mixing is very significant for violation of the CP 
symmetry in the the standard electroweak model which 
we discuss in the next chapter. The right-handed quarks 
are SU{2) singlets with hypercharge Y = 4/3 for the 
charge Q — 2/3 quarks u Ri cr and tn since the charge 
Q = —1/3 quarks d' R , s' R and b' R have hypercharge Y = 
-2/3. ' 

In a similar way the left-handed leptons are present in 
the model grouped in SU (2) doublets of the form 




where the left-handed lepton states are II = ^(1 — "f 5 )l 
and I stands for e, /i, t and the corresponding neutrinos. 
The right-handed leptons are weak isospin singlets of the 
form l R = i(l + 7 5 )Z. The I now stands for the e, /i, r, 
since it is convenient to idealize the neutrinos as massless 
and in this case right-handed neutrinos do not exist. In 
order to satisfy the Gell-Mann-Nishijima relation for the 
lepton charge the left-handed leptons have hypercharge 
Yl = —1, whilst to the right-handed we assign hyper- 
charge Y R = — 2 . 

As a conclusion to the above discussion, the fermion 
part of the Lagrangian can be written as 

+ E7fl7 M (i^-ff , ^)/fl. (13) 



where the sum runs over all left- and right-handed 
fermion fields. Addition of the gauge and fermion parts 
results in a Lagrangian density which describes massless 
gauge boson fields interacting with massless fermionic 
matter fields. 

Gauge boson mass terms which should be of the form 
\M^B^B^ plus similar terms for the A 1 ^ are clearly not 
gauge invariant. In addition a fermion mass term should 
be of the form 

-mff = -m(J R f L +J L f R ). (14) 

Such a term manifestly breaks the gauge chiral invari- 
ance, since Jl is a member of a SU(2)l doublet whilst 
f R is a singlet. Thus addition of mass terms by hand 
into the Lagrangian destroys the gauge invariance and 
the resulting theory is not renormalisable. However, it is 
an experimental fact that the vector bosons which medi- 
ate the weak interactions and the fermions are massive 
particles. The answer to the crucial problem of the mass 
generation without loosing the renormalisability of the 
theory is the Higgs mechanism, which is discussed in the 
next section. 



C. Spontaneous Symmetry Breaking 

The only known renormalisable theories involving vec- 
tor bosons are gauge theories, but since mass terms in 
the Lagrangian spoil the gauge invariance it is impos- 
sible to find a renormalisable theory of massive vector 
bosons. However, it is possible for the Lagrangian and 
the equations of motion of a system to have a symmetry 
but the solutions of the equations to violate this sym- 
metry. This phenomenon, where the solutions violate 
the symmetry of the equations, is known as spontaneous 
symmetry breaking (SSB). A well known example is the 
case of a ferromagnet, where the Hamiltonian of the sys- 
tem is rotationally invariant, but the ground state of the 
system has a non zero magnetization and it is clearly not 
invariant under rotations. The analogous of the ground 
state of the ferromagnet in particles physics is of course 
the vacuum. The Lagrangian of the theory is taken to 
be invariant under an internal symmetry, but the vac- 
uum state is not. The ground state is characterised by 
some scalar fields which are not invariant under the sym- 
metry transformations and develop a non zero vacuum 
expectation value. 

Consider the case of one scalar field with quartic self 
interaction. The Lagrangian of the theory is given by 

£ = \{d^f -ImV-^ 4 , (15) 

where m is the field mass and A its coupling constant. 
The potential energy density of this field at the classical 
level, which is called the classical potential is given by 

Vo = \m 2 <P 2 + \\<P" • (16) 
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term becomes positive 
m (cp = a) = — 



= 3Acr 2 - /i 2 = 2/i 2 > . (17) 



In determining masses and interactions one must per- 
turb about the stable vacuum by shifting the field <j) as 
4>{x) = a + f]{x) . Inserting that into the Lagrangian, 
the mass term appears with the correct sign. The above 
example serves as a useful introduction to symmetry 
restoration at high temperatures, which we discuss in 
Section EH 



D. Boson and Fermion Masses 



FIG. 1: The classical potential Vo(4>) before symmetry break- 
ing. 




FIG. 2: The classical potential Vo(<f>) after symmetry break- 
ing. 



The Lagrangian remains invariant under the the symme- 
try operation which replaces cj> by — (p. 

Let us consider now two cases and first suppose that 
the mass squared of the field is positive, m 2 = /i 2 > 0. 
From Eq. I|16|) it is clear that the potential has a mini- 
mum, which is the most favourable state of the system, 
for = 0. The minimum of the potential is shown in 
Pig- ID. 

But in the case when the mass squared is negative, 
that is m 2 = —fi 2 < 0, we observe a very different sit- 
uation . The potential energy has two minima which 
occur at <fi = a = ±/x/\/A and now there are two de- 
generate ground states. Once we choose one of them, 
the ground state does not share the symmetry with the 
Lagrangian. This situation is illustrated in Fig. |21 The 
value of the curvature of the potential near the extrema 
describes the effective mass squared of the scalar fluctua- 
tions. Prior the symmetry breaking the mass term of the 
field is negative. After breaking of symmetry the mass 



In order to understand the problem of mass generation 
for fermions and vector bosons we review here the well 
known Higgs model. In the simple scalar model, which we 
examined earlier, the mass term after symmetry break- 
ing appeared with the correct sign. In the Higgs model 
insertion of the scalar field into the Lagrangian results 
in the generation of masses for the particles which inter- 
act with. The Lagrangian of this model exhibits a U(l) 
gauge symmetry and is given by 

C = (d»+igA»)<p*(d fl -igA fl )<p-~F lll/ F>"' 

- M 2 0*0 - A(0*0) 2 . (is) 

For fi 2 > the Lagrangian, apart from the scalar self 
interaction term, is that of scalar QED. However, in the 
case where ^ 2 < 0, shifting the field to its true ground 
state as 

4(x) = -j=[a + ri(x)+it(x)] (19) 

and inserting that into the Lagrangian, the latter be- 
comes 

-X<j 2 V 2 + ^g 2 a 2 (A^ 2 ^gaA^ (20) 

where for simplicity we have omitted the interaction 
terms. The particle spectrum in the new Lagrangian 
appears to be composed by a massless scalar, the so- 
called Goldstone boson, a massive scalar particle rj with 
mass m v = \/2\a 2 and a massive vector A^ with mass 
m A = go. 

One can use the freedom of a gauge transformation 
to eliminate the unphysical Goldstone mode. Then the 
Lagrangian looks like 

C" = l -(d^h) 2 - \o 2 h 2 + l -g 2 cj 2 Al- Xah 3 -±\h 4 
+ l -g 2 Alh 2 + ag 2 A 2 t h - \f, v F^ . (21) 
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The Goldstone boson no longer appears in the theory 
and the Lagrangian describes just two interacting mas- 
sive particles. The particle spectrum consists of a massive 
vector boson and a massive scalar h, which is known 
as a Higgs particle. 

The mechanism we have presented above, is known as 
the Higgs mechanism. An overview of the model can be 
found for example in Halzen and Martin |2^| or Cheng 
and Li [3(j . The Higgs mechanism can be generalized to 
non Abelian gauge theories and in particular to the case 
of the standard model of electroweak interactions. We 
present how, due to spontaneous breaking of symmetry, 
the fermions and vector bosons acquire their masses. In 
the case of the so-called minimal standard model, one 
complex SU(2) doublet of scalar fields, with hypercharge 
Y = 1 is introduced, 



$ = 



(22) 



1 ( 4>i + i(f>2 
V2 V ^ 3 + 

where <f>, i = 1, 2, 3, 4 are real fields. The building block 
of the Lagrangian containing the scalar doublet looks like 



= (D^(D^) - , (23) 
where the covariant derivative is given by 



D, = d,~ i 9 -r a Al - i 9 -B, 



(24) 



and the most general gauge-invariant renormalisable po- 
tential is of the form 



(25) 



with the mass p? being positive. Without loss of gener- 
ality the minimum of the field $ can be chosen at 



V2 W ' 



(26) 



where a 2 = /i 2 /A and since we are interested in the physi- 
cal particle spectrum the unitary gauge is an appropriate 
choice. Just as in the case of the Higgs model, after spon- 
taneous breaking of symmetry the kinetic term of the 
scalar part of the Lagrangian (D A '$)^(Z? /i $) gives rise to 
the vector boson masses. The two charged vector bosons 
are defined by the combination 



W,, 



with a classical mass 



(27) 



(28) 



where a is the vacuum expectation value of the Higgs 
field. The neutral vector boson field appears as the com- 
bination 



9' By - gA% 



(29) 



and its classical mass is given by 



(30) 



The photon field remain massless tcia = and is given 
by the combination 



g'B^ + g A% 



(31) 



The generation of the fermion masses can be done in 
a similar way. For this we introduce into the Lagrangian 
Yukawa coupling terms between scalars and fermions of 
the general form 

C Y = g e Le^e R + g d Qu® d R + guQ u $"UR + h.c . (32) 

for all the fermion doublets. In order to generate the 
quark masses for the upper components of a quark dou- 
blet we have constructed a new Higgs doublet from <f> 



(33) 



which transforms identically to but has opposite weak 
hypercharge Y = — 1. After symmetry breaking the neu- 
trinos remain massless, while the quarks and the lower 
components of the lepton families acquire masses of the 
form 



m f 



if* 
V2 



(34) 



where gf is the Yukawa coupling of fermion / to the 
scalar doublet (f>. 

The final Lagrangian of the standard electroweak 
model appears as a sum of the gauge boson term Eq. J7J), 
the fermion term Eq. I|13|) . the above given scalar part 
Eq. 112MI) and the Yukawa interaction terms Eq. I|32l) . In 
addition to those terms a gauge fixing term must be in- 
cluded into the Lagrangian. The unitary gauge which 
we had chosen earlier is convenient and displays clearly 
the particle spectrum, but it is not the most appropriate 
when one calculates quantities involving the propagators 
of the vector bosons. Instead of the unitary gauge, a class 
of gauges like the so-called gauges can be chosen. In 
this case the gauge fixing term has the form 

C GF = - l -tgv X a ) 2 - ^s'V) 2 . 

(35) 

Different values of £ correspond to different gauges. One 
choice is the Landau gauge (£ = 0), which although has a 
bit complicated propagators, is very convenient in calcu- 
lations involving scalar potentials because the coupling of 
unphysical scalars (Goldstone bosons) to physical scalars 
is zero. On the other hand in calculations involving the 
effective potential (which we discuss in Section ITVT) the 
Landau gauge is used 0, |27]], mainly because of the 
vanishing interaction between physical scalars and ghosts 
(necessary ingredients of non abelian gauge theories). 
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E. Phase Transitions in Gauge Theories 

As it was referred in the introduction, a raising of the 
temperature results in the disappearance of the classical 
field which is responsible for the breaking of symmetry 
and this corresponds to a phase transition 0, ■ Sup- 
pose we deal with one scalar field (f>. At zero temperature, 
the location of the minimum of the potential Vq(<P) de- 
scribes the true ground state, the vacuum of the theory. 
At a finite temperature T the equilibrium state of this 
field is governed by the location of the minimum of the 
free energy density F(cf>,T) = V(<f),T), which is equal 
to Vo(<p) at zero temperature. The contribution to the 
free energy density from ultrarelativistic scalar particles 
of mass m at a temperature T which is much larger than 
the particle mass has the form 



m 2 T 2 
24 



AF = AV(4,T) = -—T 4 



O(-) 



(36) 



The resulting potential energy density at finite tempera- 
ture, which is known as the effective potential V(cf>, T), is 
the sum of the classical potential plus the above contribu- 
tion. The appearance of these contributions to the clas- 
sical potential is discussed in some detail in Section llVl 
where the theory of the effective potential is explored. As 
we have shown earlier the field dependent particle mass 
is given by 



d 2 V 



= 3A</> 2 



(37) 



so the full expression of the temperature dependent po- 
tential, at the limit of temperatures large enough com- 
pared to the particle masses looks like 



2^2 



A<?!> 4 



XT 



2^,2 



n 2 T 4 M 2 T 2 
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!)() 



24 



(38) 



At zero temperature the effective potential has a local 
maximum at (j) — and an absolute minimum at (f> = 
(7^0, but as the the temperature increases, the energy 
difference between the minimum of the potential at 4> =/= 
and the local maximum at <j) — decreases. At a critical 
temperature, which is given by 



(39) 



this energy difference disappears and the only minimum 
of V((f>, T) is the one at <j> = 0. 

This corresponds to a phase transition which takes 
place from a state with broken symmetry to a symmetric 
state and the particles appear to be massless again. The 
behaviour of the effective potential for several tempera- 
tures is given in the next qualitative picture in Fig. [3] . 

Any phase transition is described in terms of an order 
parameter which distinguishes the various phases. There 
are two types of phase transitions. In phase transitions of 
the first type the order parameter jumps discontinuously 



Vfo,T) 




FIG. 3: Evolution of the effective potential V(4>, T) in the 
case of a second order phase transition. 



from its value in the first phase to that in the second. A 
phase transition of this type is called a first order phase 
transition. Phase transitions where the order parameter 
changes continuously are known as second order phase 
transitions. In our previous example the order param- 
eter of the transition is the classical field. During the 
transition the value of the classical field at the minimum 
decreases continuously to zero and the transition is of the 
second order. 

A more complicated case appears in theories with more 
than one coupling constants, as for example the Higgs 
model, which we have discussed in the previous section. 
In the Higgs model when the temperature is high enough 
as compared to the particle masses and if one assumes 
that A ~ g 2 , temperature induced effects add a term to 
the classical potential, so it looks like 0, 0, 



V{cj>,T) 



I 

2^ 



2 j,2 



4A + 3g' 
24 



T 



2i2 



(40) 



The phase transition takes place at a critical temperature 
which is given by 



T 



12n 2 



4A + 3. 9 2 



(41) 



and we deal again with a second order phase transition, 
since the field depends on the temperature continuously. 

However, a different situation appears if we consider 
the case A < g 4 since the mass of the vector particles, 
which is given by m(</>) = g4>, can no longer be neglected 
compared to temperature and their contribution to the 
effective potential result in the form 



AV(^,T) 



3m 3 ((t))T 
12^ 



(42) 



plus field independent terms which we omit. The extrema 
of the effective potential correspond to the solutions of 
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FIG. 4: Evolution of the effective potential V(<f>, T) in the 
case of a first order phase transition. 



the equation 



- = o. 

u, v 12 4tt ' 

(43) 

The behaviour of the effective potential as a function of 
the field for several temperatures is given in Fig. 0] 

As we can observe due to the appearance of the cu- 
bic term in the field the effective potential has now three 
extrema. There are two local minima at <f> = V\ and at 
= and a local maximum at <p = v 2 ■ The state which 
corresponds to <fi = V\ is stable at low temperatures but 
it becomes metastable above the critical temperature. 
When the system reaches the critical temperature, the 
potential has two degenerate local minima, one giving 
the broken symmetry phase </> ^ and the other the 
symmetric phase cf> = 0. At the critical temperature as 
the system cools down a phase transition takes place from 
the symmetric phase to the state with broken symmetry. 
This phase transition is of the first order, since the field 
depends on the temperature discontinuously. In particu- 
lar when 3g 2 /167r 2 < A < g 4 , the critical temperature in 
the Higgs model is given by 



1/4 



(44) 



F. Cosmological Phase Transitions 

Phase transitions like the ones we presented above 
could have occurred during the expansion of the uni- 
verse. According to the standard hot universe theory, 
the universe started its life from a state with enormous 
density and temperature and with symmetry restoration 



between all the interactions. At the beginning of the ex- 
pansion all the particles appear to be massless. As the 
universe cools down phase transitions take place breaking 
the symmetry between the interactions and result in the 
generation of the massive particles we observe today. It 
is believed that the universe started from an initial state 
at a temperature which was at least of the order of the 
Planck scale about T ~ 10 19 GeV. During the expansion 
three such phase transitions could have occurred. The 
first phase transition, known as the GUT's phase transi- 
tion, could have happened at a temperature of the order 
T C1 ~ 10 14 - 10 15 GeV. Through this phase transition 
a scalar field <f> ~ 10 15 GeV is generated, which breaks 
the symmetry between the strong and electroweak in- 
teractions. As the temperature decreases further, at a 
temperature of the order of magnitude T C2 ~ 10 2 GeV, a 
second phase transition can take place and the symme- 
try between the weak and electromagnetic interactions 
breaks. This is called the electroweak phase transition 
and is the main subject of this work. The nature of the 
electroweak phase transition will be investigated in detail 
in Section Finally, at a temperature T ~ 10 2 MeV, 
the QCD phase transition takes place which breaks the 
chiral invariance of the strong interactions. 



III. BARYON NUMBER AND CP VIOLATION 



Introduction 



As it was stated in Section [i] the Sakharov conditions 
demand that baryon number and CP must be violated, 
in order to produce a net baryon number. Since the work 
of t'Hooft |23, baryon number violation has been shown 
to occur in the standard electroweak model through the 
axial anomaly, although the rate of anomalous baryon 
number non conserving processes are exponentially sup- 
pressed at zero temperatures. As the temperature in- 
creases the rate of baryon number violation increases and 
is no more negligible at temperatures close to the elec- 
troweak scale. The electroweak theory violates charge 
conjugation C, while charge conjugation combined with 
parity (CP) is violated through the quark Yukawa cou- 
plings, although this CP violation appears to be too 
small to generate the required asymmetry. 

Our aim is to sketch how baryon number violation and 
CP violating effects take place within the physics of the 
standard electroweak model. Baryon number violation at 
zero and finite temperature is discussed in Section IlII Bl 
while the final Section IIII Gl is devoted to CP violation 
as it appears in the standard electroweak model. 
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B. Baryon Number Violation 

1. Zero Temperature 

In the standard electroweak model baryon number J% 
and lepton number j£ currents are exactly conserved 
at classical level. However quantization of the theory 
leads to the appearance of anomalous axial currents [2(| • 
In what follows we ignore electromagnetic effects, which 
mathematically correspond to setting the Weinberg angle 
to zero 0, HE El ■ It is as if by quantum corrections 
the W and Z boson fields have acquired an indeterminate 
baryon plus lepton number. 

Violation of the baryon number and lepton number 
currents in the standard electroweak model comes from 
two crucial facts. The first is that due to the nature of 
the electroweak interactions the gauge boson fields couple 
differently to the left-handed and right-handed fermion 
fields. When a fermion couples to a gauge boson the 
corresponding axial current is anomalous 

d^ = N F ^F^F a ^ , (45) 

where Np is the number of fermion generations and g is 
the SU(2)l coupling constant. The field strength tensor 
of the SU(2)l group is given by 

F; v = dpAl - d v A« - ge abc A;al , (46) 

and is related to F^ v through 

F„ v = \e^ p °F pa . (47) 

In the above equation F pa = F pa T a where T a are the 
generators of the gauge group algebra. Calculation of the 
anomaly for the lepton number current gives an identical 
result, so B — L is conserved but B + L is not. 

It is very significant that the term on the right hand 
side of the Eq. may be written as four-divergence of 
a new current, the so-called Chcrn-Simons current 

JC" = ^ vp °(F* p A a a - -e abc A a v A b p A%) (48) 
3 

since one can find that 

d M K" = F^F a " u . (49) 

The second crucial fact is that the standard elec- 
troweak model with an SU (2) gauge group and a Higgs 
field possesses a non-trivial vacuum structure. In addi- 
tion to the usual vacuum configuration A{ = A p T a = 
and <fi = 4>o, there is an infinite number of pure gauge 
transformations of the form 

4>' = r(n)to 

A'- = -diSlQ- 1 , (50) 
9 




1/2 3/2 5/2 7/2 

H, 



FIG. 5: The potential energy of the SU(2) gauge field as 
a function of the Chern-Simons number. The minima cor- 
respond to configurations with zero gauge field energy but 
different baryon number 



where Cl(x) is a 2 x 2 unitary matrix and T(f2) an ar- 
bitrary rotation matrix and under these transformations 
the Lagrangian of the model remains invariant. 

Physics imposes constraints, so one can find a non-local 
variable but locally gauge invariant quantity defined as 

Ncs = J^e yfe J d 3 x(F%A a k - 2 -e abc A1A)Al) , (51) 

which characterizes these degenerate vacua and it known 
as the Chern-Simons number. Under "small" gauge 
transformations the integrand of Nqs changes by a quan- 
tity which is a total divergence. But under a "large" 
gauge transformation £1 — > 1, as the fields approach zero 
at the spatial infinity, the Chern-Simons number Ncs is 
shifted by an integer number n which is called the wind- 
ing number. The winding number is given by 

which is a non-trivial topological integral and it can be 
any real number but for vacuum configurations of the 
field it is an integer. The winding number characterizes 
the degree of mapping of the space compactified as a 3- 
sphere onto the group SU(2). 

The classical potential energy E of the SU(2) gauge 
field as a function of the Chern-Simons number Ncs is 
multiply periodic and it can be schematically represented 
as in Fig. |3J 

Between the degenerate vacua there exist unstable, 
time-independent solutions of the field equations called 
sphalerons. The sphalerons were the first introduced 
by Klinkhamer and Manton and they have Chern- 
Simons numbers Ncs = 1/2,3/2.... At these saddle 
points the height of the barrier separating topologically 
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distinct vacua has been estimated in the range 

E sph = ^L B {\/g 2 ) =8-14 GeV (53) 

for A varying from zero to infinity. 

It is very significant that this non-trivial vacuum 
topology is related to the anomaly equation. The total 
baryon number is given by 

B = J d 3 xJ% (54) 

and it can be shown that the baryon number change 
during any time interval is related to the Chern-Simons 
number through the equation 

AB = N F [Ncs(h)-N C s(ti)} , (55) 

since the time component of JsT M in Eq. 148|) is simply 
the Chern-Simons number Ncs- Therefore if the baryon 
number is not to be conserved it is related to a change in 
the Chern-Simons number of the SU(2) gauge vacuum. 

At zero temperatures or low densities a transition from 
one vacuum configuration to another can be achieved by 
quantum tunneling. In the computation of the transi- 
tion rates, a classical solution of the equations of motion 
can be used which is called an instanton. The instantons 
correspond to the classical solutions used in the WKB 
method of quantum mechanics |3lL 15(1 |5jJ. The transi- 
tion rates have been computed by t'Hooft[l^.l20|. but as 
it was shown these are suppressed by a semiclassical fac- 
tor exp(— Air/aw) ~ 10~ 150 . At zero temperature there- 
fore baryon number non conservation is negligible. This 
of course is a consequence of the very large energetic bar- 
rier (~ 10 GeV) separating the vacua of different baryon 
number due to the sphaleron but fortunately this picture 
is completely different at high temperatures. 

2. Nonzero Temperature 

At finite temperatures thermal fluctuations with an en- 
ergy higher than the height of the barrier between the two 
degenerate minima can classically traverse the maximum 
and enter the new minimum, leading to unsuppressed 
baryon number violation. The rate of penetration of the 
sphaleron barrier is then given by the Boltzmann factor 
associated with the formation of the sphaleron 

r = (%T) 4 eM-E sp h{T)/T] , (56) 

where the sphaleron energy at temperature T is given by 

E sph ~ B^l . (57) 
aw 

The factor B is obtained numerically by solving the field 
equations. Such a numerical solution is still lacking, how- 
ever we can expect that for temperatures larger than the 



W boson mass, the baryon violation rate should behave 
like 

T - awT 4 , (58) 

due to the vanishing of masses at high temperatures. 
This of course leads to a completely unsuppressed rate 
of baryon number violation. In particular this means 
that any pre-existing baryon asymmetry generated dur- 
ing the GUTs phase transition, would have disappeared 
until the time when the EW phase transition takes place. 
However these processes only violate B + L since as we 
have seen B — L does not have an anomaly and is exactly 
conserved in the standard model. Therefore a B — L ex- 
cess produced in the early universe will not be washed 
away. 

It is very significant the role played by the sphaleron 
after the completion of the electroweak phase transition 
[33l l36l |. The sphaleron mass, which is proportional 
to the Higgs vacuum expectation value a(T c ), must be 
large enough to strongly suppress baryon number vio- 
lating processes. Otherwise, these processes may erase 
any asymmetry produced during the phase transition. 
For baryons to survive we need the sphaleron rate to be 
much less than the expansion rate of the universe. These 
considerations impose the constraint 

a(T c )/T c > 1 , (59) 

where T c is the temperature immediately after the tran- 
sition and can also be used in order one to obtain an 
upper bound for the Higg s mass, which has found to be 
m H < 45 GeV jiJEiif. 

C. CP Violation 

It has been well understood that C and CP violation 
is a crucial feature of any theory that attempts to explain 
the observed asymmetry between matter and antimatter 
in the universe starting from initially symmetric condi- 
tions. CP symmetry requires that the partial rates of 
CP-conjugate processes be equal, therefore for any pro- 
cess that violates the baryon number there would be a 
CP-conjugate process of equal rate and no asymmetry 
could be generated. So far the Kaon system has been the 
only laboratory for the observation of CP violation. 

Violation of C symmetry is a generic feature in the 
standard model of electroweak interactions, since due to 
the chiral nature of the SU (2) gauge interaction, left and 
right chiralities of quarks have different interactions, so 
the Lagrangian of standard model has both C-even and 
C-odd pieces. 

The standard model with two fermion families and a 
single Higgs doublet is a CP invariant theory, but addi- 
tion of the third generation violates the CP invariancc 
of the Lagrangian. As we referred in the second chapter 
on the formulation of the minimal standard model with 
one Higgs doublet, CP violation in the model originates 
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from the Cabbibo-Kobayashi-Maskawa (CKM) matrix 
which relates the quark mass eigenstates to the weak 
eigenstates. 

The Yukawa term of the standard model Lagrangian 
involving the quark interactions with the scalar doublet 
was given in Eq. (|32|l . but in terms of fields which are 
eigenstates of the weak interactions, can be written in 
the form 



9w 



V2M W 



Q L KM d d a $Q L KM u u R $ + h.c. 



(60) 

where K is the CKM matrix which relates the two set of 
eigenstates and M u , M d are the diagonal mass matrices 
of the quarks with charges 2/3 and —1/3 respectively. 
The specific form of the CKM matrix can be found in 

For three fermion generations the CKM matrix must 
be 3 x 3 and unitary and this constraint provides relation- 
ships between its elements, so finally its parameterisation 
can be done through three rotation angles and there can 
be only a physically meaningful phase which leads to ob- 
servable CP violation effects HjEIH- This CP vio- 
lation can be quantified through the combination |45| 



A 



CP 



sin t/12 sin #23 sin W13 sin ocp 

,2 „,2w 2 rm 2\( rm 2 2 



sin 5c 

x(mf - m z c )(mf - m z u )(m 2 c - m z u ) 
x(ml-ml)(ml-ml)(m 2 a -m 2 d ) , (61) 



where fty, i, j = 1, 2, 3 are real angles and the phase Sep 
lies in the range < 5c p < 2-7T. In order to get nonzero 
CP violation, the quark mass matrices must have non- 
degenerate elements and also 9ij ^ 0, n, i,j = 1,2,3, 
while 5c p =/= 0, tt. 

According to the analysis by Shaposhnikov [42|. the 
strength of the CP violation in the standard electroweak 
model even at high temperatures appears to be too small, 
so it is unlikely that it can explain the observed baryon 
asymmetry. This conclusion results from arguing that 
the only natural scale for the baryogenesis problem is the 
temperature of the electroweak phase transition, which as 
we show in the Section^is of an order T ~ 100 GeV. At 
this temperature the Yukawa interaction can be treated 
as a perturbation, because the quark masses are small 
compared with the temperature. The baryon asymmetry 
as defined in Section [fl is a dimensionless number, so to 
get an estimation of the asymmetry, the quantity given by 
Eq. H61|) should be divided by something with dimensions 
of (mass) 12 and a natural choice is the temperature itself 
H3. Hence, 



n B 



< 



A 



CP 



NeffT 



12 



10 



-20 



(62) 



where N e ff is the number of degrees of freedom in the 
standard model. 

The same author suggested in subsequent work 0,0] 
a possible way out of this problem. He assumed that 
accumulation of small CP violating effects during the 



expansion of the universe could enhance the total CP vi- 
olation during the electroweak phase transition, but this 
result is not acceptable in general. Recent investigations 
show that the most promising models are models with 
enhanced the scalar sector and in particular the models 
with two scalar Higgs doublets |5flj . 



IV. THE EFFECTIVE POTENTIAL 
A. Introduction 

At the very early stages of its evolution the universe is 
filled with matter at very high energies and densities, so 
matter should be described in terms of quantum fields. In 
order to investigate how the electroweak phase transition 
could have occurred at the beginning of the universe we 
need to work in the framework of quantum field theory. 
The basic tool for the investigation of the nature of the 
electroweak phase transition is the effective potential, the 
quantity which has the meaning of the potential energy 
density of the system under consideration. 

The discussion on spontaneous symmetry breaking in 
Section Til CI was purely classical. The particle spectrum 
was determined by minimization of the classical potential 
Vo(<j>) as it appears into the Lagrangian and describes the 
potential energy density of a constant scalar field. The 
effective potential has also the meaning of a potential 
energy density. A quantum field theory involves virtual 
particles, which affect the field energy density through 
emission and reabsorption processes. This generalization 
of the classical potential to include the quantum correc- 
tions is known as the effective potential. Minimization of 
the effective potential gives the field configuration with 
the minimal energy, the vacuum of the theory. 

Proceeding further, analysis of matter behaviour at 
non zero temperatures involves thermal fluctuations of 
the fields that one should take into account. Thus, a gen- 
eralization of the effective potential at finite temperature 
is needed, for the inclusion of temperature dependent 
quantum effects. As it will be clear in what follows from 
the mathematical definition the effective action has the 
meaning of the free energy of the quantum system under 
consideration. The finite temperature effective potential 
V(4>, T) as Linde states H3 at its extrema coincides with 
the free energy density. 

The effective potential has been studied extensively 
in the literature. An elegant discussion on the physi- 
cal meaning of the effective potential and its calculation 
is explored by Coleman [5(j ■ A detailed analysis of the 
theory of the effective potential at zero and finite temper- 
ature with applications on cosmological models is given 
by Brandenberger [fS^] ■ The electroweak Higgs potential 
for the standard model and its extensions has been in- 
vestigated by Sher [27}]. Early discussions on the subject 
are those of Coleman and Weinberg [4!j involving calcu- 
lation of the effective potential for a general gauge theory 
at zero temperature. Generalization of the effective po- 
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tential at finite tem per ature is given by Dolan and Jackiw 
[l7| and Linde 0, [2^ . In what follows we give a formal 
discussion on the notion of the effective potential as it 
appears in the framework of quantum field theory. 

The outline of this section is as follows: In Section llV Bl 
we discuss the properties of the effective potential at zero 
temperature and give the basic steps for its calculation in 
the framework of quantum field theory, using the path in- 
tegral formalism. For this analysis we restrict ourselves 
to the one loop radiative corrections. In Section II V CI 
we generalize the idea of the effective potential at finite 
temperature. This generalization appears in a form of 
integral terms which we add to the zero temperature ef- 
fective potential. At high temperatures these integrals 
can be approximated expanding them in a series. The 
high temperature approximation of the fermion and bo- 
son contribution to the effective potential is compared 
with the exact calculation of the one loop radiative cor- 
rection integrals using numerical methods. 



B. Zero Temperature Effects 

1. Path Integral Formalism 

The effective potential can be calculated in an elegant 
way using the path integral formalism and the notion of 
generating functionals. Consider the case of a real scalar 
field and suppose that an external c-number source J(x), 
a function of space and time, is added into the Lagrangian 
coupled linearly to the field <f>. Then the transition ampli- 
tude from the vacuum state in the far past to the vacuum 
state in the far future is defined as 



V(j)exp(i / d A x[C{(j)) + J(x)<p(x)]) . (63) 



Z[J] 



This amplitude can be expanded in a functional Tay- 
lor series in terms of the source J(x), with coefficients 
G™(xi, X2, ■ ■ ■ , x n ) the Green functions of the theory. 
These Green functions can be found by functional dif- 
ferentiation of Z[J] with respect to J(x), at J(x) = 0. 

One can introduce the generating functional of the con- 
nected Green functions W[J], which is related to Z[J] by 



Z[J] = exp(iW[J] 



(64) 



Expanding this functional in powers of J(x), the coeffi- 
cients G™(xi, X2, ■ ■ ■ , x n ) are the connected Green's func- 

(n) 

tions. The connected Green's functions G c are given by 
the functional derivative of W[J] with respect to J(x), 
at J(x) = 0. 

At this point we have to introduce a quantity which is 
referred as the classical field <fi c (x). This is a functional of 
the source J(x) and is defined as the vacuum expectation 
value of the field operator in the presence of the source 
as 

SW[J] 



4>c{x) 



SJ(x) 



(65) 



In the absence of the source, the classical field <p c (x) is 
just the vacuum expectation value of the field operator. 

The generating functional of the 1PI (one particle ir- 
reducible) Green functions is called the effective action 
and is defined by the functional Legendre transform 



r[0 c ] = w[j] - / d 4 xJ(x)Mx) 



(66) 



It is called the effective action, because it is a functional 
of the classical field 4> c and hence akin to the classical 
action S [(/>]. From the definition of the effective action, 
it follows that the source can be obtained in the form of 
a functional derivative as 



54> c (x) 



(67) 



The effective action can be expanded in powers of <j) c , or 
alternatively one can expand about a constant value of 
the field <f> c . This is the same as expanding in powers 
of the derivatives of <f> c . Thus in position space it is an 
expansion of the type 



1 



d*x[-v(<j) c ) + -(d^ c yz(t c ) + 



(68) 



where V(<p c ) and Z(<p c ) are ordinary functions of <f> c , not 
functionals, and the function V(<p c ) is known as the ef- 
fective potential. In the case of a classical field <j) c which 
is constant in space and time and in the absence of the 
source J, <p c has the significance of the vacuum expecta- 
tion value (VEV) of the field operator 



dV(</> c ) 
d(p c 







(69) 



This last equation allows one to interpret the VEV of 
the field operator as the stationary point of the effective 
potential, which can be obtained by solving the above 
equation. Of course, in the case when the field is to 
vary in space or time, we have to solve the more general 
equation ST/S(j) c = 0. Thus the VEV of the field op- 
erator, taking into account quantum corrections, can be 
found by minimizing the effective potential. If the effec- 
tive potential has several local minima, it is the absolute 
minimum that corresponds to the true ground state. 

It is useful to see that the effective potential can be 
obtained as the infinite sum of the 1PI graphs with van- 
ishing external momenta |49| , although we will not make 
use of this method in order to evaluate that in the case 
of the scalar field. One can expand the effective action 
functionally in terms of 4> c as 

°° 1 f 

r kM = ^ — / d 4 xi . . .d 4 x n 4> c {xi) . . . <j> c (x n ) 

n=l n ' 

xT^( Xl ,...,x n ) . (70) 



The coefficients are referred as the one-particle- 

irreducible (1PI) Green functions. Fourier transforming 
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the 1PI Green functions to momentum space, the expres- 
sion for the effective action takes the form 



r[0 c ] = 



1 

n=l 



d 4 pi . . . d 4 p n 5 4 {pi + . . . + p n ) 



xT^{p 1 ,...,p n )^ c (p 1 )...M Pn ) , (71) 



where <p c (jPi) is the Fourier transform of <p c . Thus, com- 
paring this last expression with the momentum space ex- 
pansion of the effective action above, we find that the 
effective potential is given as a sum of the 1PI graphs of 
the form 



-rW(ft = o) <p n c 

n=l 



(72) 



2. Scalar Loops 

In order to understand how the effective potential is 
calculated, we present here the evaluation of the effective 
potential for a theory involving a real scalar field with 
quartic self interaction. The Lagrangian governing this 
theory is given by 



(73) 



The potential energy density of the real scalar field at 
classical level, which is commonly called the classical po- 
tential and denoted by Vq(<W, is given by 



F o (0) = -^V + 



(74) 



where the field mass squared /i 2 and the coupling con- 
stant A are positive. The effective potential which is usu- 
ally denoted by V(<fi) will result in as a sum of the clas- 
sical potential, plus terms due to radiative corrections 

As it was referred earlier the diagrammatic method of 
summing the 1PI graphs can be used for the calculation 
of the effective potential, however we present here an 
alternative method which is based on the saddle point 
evaluation of path integrals 25, 26] . According to this 
method exponential integrals involving a function f{x) 
which is stationary at some point xq, can be solved by 
expanding f(x) about this point. Then, omitting higher 
derivatives, the integral becomes a Gaussian and thus 
can be evaluated according to standard methods. 

The starting point is the generating functional of the 
connected Green's functions, where we have restored the 
Planck's constant so it takes the form 



exp ^W[J]j =Af J 2?0exp (±S\ 



(75) 



and the normalization constant J\f is chosen such as 
W[0] = 1. Although in natural units the Planck con- 
stant is equal to unity, we restore it here in order to get 



an expansion in terms of h and make clear that we calcu- 
late quantum corrections to the classical potential. The 
powers of H count the number of the closed loops in the 
loop expansion. The action S[cj>, J] in presence of the 
source J(x) is given by 



S[<f>, J] = J d 4 x[£(0(x))+J{x)(f){x)} 

The saddle point is at <f> — 0o an d satisfies 
SS[i>, J] 



5<f>(a 



h,j(x) 



(76) 



(77) 



where cj> is a function of x and also a functional of J. In 
the case when the source J tends to zero, 4>o(x) becomes a 
solution to the classical equations of motion. Expanding 
the action S[<j>, J) about the stationary point at 4> = (f>o 
yields 



S[4>,J\ = S[^o,J]+ J dx[^{x)-US'[h,J] 

dxdy[(j)(x) - 0o][0(y) - 0o]<5"[0 o , J] , 



plus higher order terms which we omit, since we are inter- 
ested in calculating radiative corrections up to one loop 
only. In this last equation we have used the shorthand 



ss 



S"[<Po, J] 



5(j>{x) 
6 2 S 



8<f>(x)5<f>(y) 



4>=<t>o 



4>=4>o 



We then differentiate functionally the action and insert 
the result into the equation above. The resulting integral 
is a Gaussian one, so going over to Euclidean space we get 
finally the loop expansion for W[J], ignoring correction 
terms of order h 2 and higher, 



w[j] = s[<i>o] + h 



dxMx)J(x) + y Tr ln[D + Vtf'fao)] 



. (78) 

Inserting this expression into the defining equation of 
the effective action we get, by putting the source J — > 0, 



I> c ] = S[4> c ] + yTrlnp + <(0 O )] • (79) 

For a field which is constant in space and time, through 
the expansion of the effective action we find for the effec- 
tive potential the expression 

h f dk 4 

V(4>) = Vo(cb) + - J ln[fc| + Vi'M] , (80) 

where in the final step we have used that the trace of an 
operator is the sum over its eigenvalues and we expressed 
the result on going over to Euclidean momentum space. 
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3. Renormalization 

The above integral is divergent as it happens in general 
when one calculates radiative corrections involving inte- 
grations over the internal momenta of the graphs. For 
the theory to be renormalized, the divergences must be 
absorbed, if possible, into the parameters of the theory. 
In order to evaluate this integral we introduce a cut-off 
at some large momentum k 2 — A 2 , so we obtain 



= Vote) 



A 2 



-m 2 {(t>) 



32ir 2 

™ 4 W>) I", (m 2 {(t)) 



64tt 2 



hi 



A 2 



Vcttt) ,(81) 



where the field dependent squared mass of the scalar 
m(<f) is defined as the second derivative of the classical 
potential as 



m 2 (</>) 



d 2 V (cf>) 



= 3\<j> 2 - fjj 



(82) 



To remove the cut-off dependence we have introduced a 
counterterm potential which has the same structure as 
the original potential 



1 



KtW>) = -A<jy z + -B^ + C 



(83) 



We can determine the coefficients A and B by requiring 
that the position of the minimum of the effective poten- 
tial and the Higgs field mass remain in their classical 
values [l(J,|23], so 







(84) 



and the Higgs field mass m 2 H = 2/j 2 results as the second 
derivative of the potential evaluated at <f> c = a, 



dV 2 ^) 



2n 2 



(85) 



The final expression of the effective potential for the 
real scalar field, including radiative corrections up to one 
loop, takes the form 



KM 
1 



TO 4 (</>) 

64tt 2 



In 



2 (0) 



m 2 (a) 



32tt 



m 2 (<f>)m 2 (a) 



128vr 



rm 4 W , (86) 



where Vote) is the classical potential. Other renormal- 
ization conditions can also be introduced [2?l |49| but the 
physical results are meant to be the same. 



effective potential in the case of the one scalar field. In 
order to proceed and obtain an expression of the effective 
potential for the standard electroweak model we need to 
know the one loop contributions of fermions and vector 
bosons. The case of a general non-abelian gauge theory 
was the first discussed by Coleman and Weinberg |49| . 
An analytic presentation of the subject is also given by 
Rivers |22| and Sher |27j. In this general case it is more 
convenient to revert to the definition of the effective po- 
tential from the expansion of the effective action and the 
first step is to extend the scalar sector including more 
fields. A different approach, based on the evaluation of 
Gaussian integrals as in the case of the scalar field, is 
given in Bailin and Love |23j |. 

According to Anderson and Hall 53], the analysis 
given in the above references can be summarized into 
that the unrenormalized one loop self-energy contribu- 
tion of virtual particles, adds to the classical potential a 
term which in the case of fermions is given by 



+ 



128tt 2 fKV > 



m 2 ((/>) 
A 2 



32tt 



2 ™}(<t>W 



(87) 



for each fermionic degree of freedom. On the other hand 
the bosons contribution has a similar form, except the 
opposite sign, and is given by 



1 



64tt 2 



1 



1287T 2 



mite) ln 

\te) 



m 2 b ((f>) 
A 2 
1 

" 32^ 



■mli^A 2 . (88) 



In both the above equations A is a cut-off which has been 
introduced in order to evaluate the divergent integrals. 

Summarizing the above discussion, the one loop effec- 
tive potential in the case of a general gauge theory in- 
volving fermion, scalar and boson fields can be written 
in the general form 



w) = vo{<i>) + wte) + vM 



(89) 



where AVi(0) is the sum of fermions, scalars and bosons 
one loop contributions to the effective potential and 
Vctte) i s a counterterm potential which is used to re- 
move the cut-off dependence as in the scalar field case. 
We adopt the same renormalization condition as before 
in order to retain the Higgs mass and the position of the 
minimum of the potential for each degree of freedom to 
which the scalar couples, so for the bosons we have a 
contribution of the form 



4- Fermion and Boson Loops 

The saddle point evaluation of the one loop corrections 
to the effective potential was an elegant way to obtain the 



AV lb = ]T 



64tt 2 



|r»in 



+2m 2 b te)m 2 (a) 



Ite) 



(°0 



\te) 



(90) 
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while for the fermions the contribution reads as 
AV lf 



E 



64tt 2 



mf(4>) In 



3 

(91) 



+2m}{cj>)m}{(j) 



where rib and n / are the numbers of degrees of freedom 
associated to spin, particle-antiparticle states and inter- 
nal symmetries (coloured quark states). 



C. Finite Temperature Effects 

1. Fields at Finite Temperature 

In order to investigate the electroweak phase transition 
in the next section, we need a generalization of the effec- 
tive potential at finite temperatures. We present here 
how the notion of the effective potential can be general- 
ized at finite temperature, following the ideas of Dolan 
and Jackiw H3, Linde [TUlHI and Sher [27|. 

At finite temperature a field theory is equivalent to 
an ensemble of finite temperature Green functions. The 
average value of an operator at nonzero temperature is 
defined by the Gibbs average as 



0(x%,x 2 , ...,x n ) 



Tre-P H 0{ Xl , 



x 2 



(92) 



where H is the Hamiltonian of the system under consid- 
eration and (3 = l/kgT = 1/T, since for the sake of sim- 
plicity the Boltzmann's constant kg can be taken equal 
to unity. Then, according to the proof given by Sher 
|27| , Green's functions at non zero temperature obey the 
same equations as those at zero temperature, but un- 
der different boundary conditions. The finite temper- 
ature Green's functions concerning bosons are periodic 
in Euclidean time, r = it, with a period (3 — 1/T, in- 
stead of having the usual boundary conditions t = ±oo. 
On the other hand fermionic finite temperature Green's 
functions obey antiperiodic boundary conditions with the 
same period (3. 

In analogy with the previous section the finite tem- 
perature effective potential can be calculated by using 
similar methods as in zero temperature field theory and 
the path integral formalism. A finite temperature effec- 
tive action r^[<^ c ] is defined by analogy to that at zero 
temperature, since the vacuum expectation value of the 
classical field at zero temperature now corresponds to a 
thermodynamic average. The finite temperature effective 
potential may be defined by an expansion of the effective 
action analogous to that of the previous section and its 
calculation proceeds through the evaluation of Gaussian 
path integrals at Euclidean space-time. The only dif- 
ference in these calculations is that when one calculates 
integrals involving boson fields all the boson momenta 
ko should be replaced by the Matsubara frequencies for 



bosons 2imT with n an integer. On the other hand, 
since fermionic fields are antiperiodic, the corresponding 
fermion momentum ko must be replaced by (2n + l)irT. 
So finally, instead of integrating over kg , one has to sum 
over n. We show this calculation in some detail and give 
relevant references in the following sections. 



2. Scalar Fields at Finite Temperature 

Consider first the case of a scalar field. We have to 
expand the Lagrangian around a constant field 4> c . The 
field dependent mass squared, which is called the effective 
mass, is given by 



m 2 (</>c) = 3A^ _ M 2 



(93) 



The zero loop effective potential, the so-called classical 
potential, is temperature independent and is given by 



(94) 



We have already calculated the one loop approximation 
to the effective potential at zero temperature as 

VfM = \J Mfc 2 +m 2 (0 c )] . (95) 

At finite temperature the above expression, as we have 
stated earlier, one has to replace the scalar boson mo- 
menta kg by the Matsubara frequency 2-KnT ', so it be- 
comes 

V?Mc) = \T £ /^ln[k 2 + (2™T) 2 + m 2 (0 c )]. 

n— — 00 

(96) 

The sum over n diverges, so in order to evaluate this in- 
tegral one can follow the procedure given by Dolan and 
Jackiw |17|. The result splits in two parts, one tempera- 
ture independent part 



d 3 k \/k" + m 2 (0 c ) 
(2tt) 3 



(97) 



which is equivalent to the one loop effective potential at 
zero temperature which we have calculated in the previ- 
ous section, and a temperature dependent part 



2^ 



dx x 2 ln[l — exp(— A(x, m, T))] 



In this last equation we have introduced the shorthand 
expression A(x,m,T) — y/x 2 + m 2 (</> c )/T 2 . These final 
steps are explained in some detail in the next section 
where we give the basic steps for the calculation of the 
finite temperature effective potential for a general gauge 
theory. 
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Thus summarizing the above discussion, the finite tem- 
perature effective potential for a real scalar field including 
the one loop radiative corrections can be written as 



V ls (4>,T) = V {4>) + Vu(4>) + v£(4>,T) 



(99) 



The finite temperature contribution vanishes as it should 
at zero temperature T — 0, when the mass squared 
m 2 (4> c ) is p ositive. But as it was stressed by Dolan and 
Jackiw if one is to use this full expression of the 
one loop potential, there are some serious difficulties for 
small values of the field, 4> 2 < /i 2 /3A, since the mass 
squared of the Higgs scalar becomes negative leading to 
an unacceptable complex one loop effective potential. 

Fortunately in the case of the minimal standard model, 
which we discuss in the next section, the scalar loops can 
be safely ignored, since they are negligible as compared to 
fermion and vector boson contributions, simply because 
there is a large degeneracy factor associated to the latter, 
while there is only one Higgs scalar. 



appears as 



TrlnC = 2 / d*x £ / ^^^+(41 ■ 

r ^ ' n 

(103) 

These can be evaluated by using the Matsubara fre- 
quency sums which in the case of bosons is given by 



Eh 



Lol +X 2 ) 



Six + 21n(l 



-0X\ 



(104) 



plus an ^-independent constant. In the fermion case the 
relevant formula has a similar form apart from a minus 
sign, 



£ ln (<4 + x 2 ) = [3x + 2 ln(l + e-P x ) 



(105) 



where for the antiperiodic fermions lu h — (2n + l)7rT. 
After summing, the scalar contribution takes the form 



3. Fermions-Bosons at Finite Temperature 

The above discussion can be applied to a general gauge 
theory and as we stated earlier there is an analogy be- 
tween the computation of the effective potential at finite 
temperature to that at zero temperature. We have to 
expand the scalar fields around their expectation values, 
which are now thermal averages, and isolate the terms 
in the Lagrangian which are quadratic in all the fields. 
Then the one loop contribution to the effective action is 
obtained as a Gaussian path integral. Suppose we deal 
with a general gauge theory involving fermions, scalars, 
gauge bosons and Fadeev-Popov ghost fields. We in- 
sert the quadratic Lagrangian into the expression of the 
effective action and evaluating the resulting Gaussian in- 
tegrals, the finite temperature effective potential results 
as a sum of traces of the form [23| 

- f dr [ (fxVftec) = -iTrlnA--TrlnB+TrljiC. 
Jo J 2 2 

(100) 

The scalar contribution is given by the first trace 

„ scalar „ 73, 

TrlnA = J d 3 X £ J j^s £ H"l + ^ + . 

(101) 

where m 2 = m 2 ((j) c ) are the eigenvalues of the scalar 
mass matrix. The contribution of the vector bosons has 
the form 

r boson r 

TrlnB = J d*x J ^ £ 3 ln[^ + k 2 + (m 2 ) a ] , 

(102) 

with m 2 = m 2 ((f> c ) being the eigenvalues of the vector 
boson mass matrix. In a similar way the fermion term 



AV ls 



a) 2 J (2n f^ 



In [1 - exp[-/3(E 2 ) i ]] 



(106) 



the vector bosons contribute the following term 

d 3 k 



AV lb (0 c ,T) 



1 r d^k v ^ 

2 J (2tt)3 2- 



6 



+ -ln [1 - exp[-/3(£; 2 ) a ]] 



and the contribution of the fermions reads as 

d 3 k 



(107) 



AV] 



if {9c 



T) 



(2ny 



£ 



(E 2 f)r 



+ -ln [1 + exp{-f3(E 2 f ) r ] 



(108) 



where we have used that E 2 = yk + m 2 for a = s,b,f 
for scalars, boson and fermions. All the above equations 
can be separated into a part AV^(0), i = s,b, f which is 
temperature independent and a temperature dependent 
part AV^(<fi) for each of the above field contributions. 
By using the fact that [H 



d 3 k 

(2^)3 



Vk 2 + R 



d 4 k 



ln(fc 2 + k 2 + R) (109) 



plus a constant independent of R and evaluating the in- 
tegrals in the zero temperature part, we get the one loop 
quantum corrections which we have already calculated 
in the previous section. By setting x = \k\/T in the 
temperature dependent part of the above equations, the 
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scalar field contribution to the effective potential at finite 
temperature takes the form 



2^2 



dx x 



x V in [i - ex P (- V^TRyr^)] , 



the boson term is 



AV^T) = i— J dxx 2 

x In [l — exp(- 



V^+K)^)] . 



and the fermion term is 
AV^T) = -4— j / rf« 2 

x ^ In [1 + exp(- ^ + (m 2 ) r /T 2 )] , 

r 

where = mi(0 c ), i = s, f,b are the scalar, fermion and 
boson mass eigenvalues. 

It is easy to observe the similarity in the formulas be- 
tween the vector bosons and scalars. The only difference 
is the factor 3 which appears in the boson case. This fac- 
tor expresses the one longitudinal and the two transverse 
degrees of freedom of a massive vector boson. On the 
other hand the factor 4 in the formula for fermions cor- 
responds to the two fermionic spin states times the two 
particle antiparticle states. On summarizing the above 
discussion we can say that finite temperature effects for 
bosons takes the form 



AV£(*,T) 



x ln[l - cx P (-^/.x 2 +to 2 (0)/T 2 )] 



The fermion term has a similar form apart from the posi- 
tive sign in the argument of the logarithm and that their 
contribution has an overall minus sign, so 

AV-(4>,T) = -J2^fdxx* 

x ln[l + exp{-^x 2 +m 2 (»/T 2 )] . 

As it was referred above for massive vector bosons the 
relevant factor which corresponds to the particle helicity 
states is = 3, since for fermions nf = 4. For coloured 
fermions, as it is the case of quarks, nf = 12. 

4- High Temperature Approximation 

It is sometimes convenient to approximate the above 
integrals by using the high temperature expansion. When 
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FIG. 6: A fermion's contribution to the effective potential as 
a function of mf(<f>)/T. The solid line represents the exact 
calculation, while the dashed line corresponds to the high 
temperature expansion. 



the temperature is high enough as compared to the par- 
ticle masses the above integrals can be approximated by 
expanding them in a series in powers of m{<j))/T. This 
expansion has been extensively used in the literature in 
investigations of the effective potential. The expansion 
of the integrals in the high temperature limit can be done 
by using the Riemann Zeta function [TtI 0] ■ According 
to the analysis given in the above references, the integral 
in the fermion case can be expanded as 



AV X H S T {4>,T) 



/ 

n 

+ 



7tt 2 T 4 m 2 f T 2 
720 + 48 



64tt 2 



In 



2^2 



O 



\T 2 



where the constant Cf is defined as In c/ = | + 2 In ir — 
27 w 2.64 and the Euler constant is 7 w 0.577. The 
particle masses which appear in the above expression are 
field dependent, so nif = m/ (</>). 

Our aim was to verify the reliability of the high tem- 
perature expansion, so we have calculated the integrals 
given in the previous section numerically. For this calcu- 
lation we have used a numerical code based on Simpson's 
rule. 

In order to compare a fermion's contribution to the ef- 
fective potential by using the high temperature approx- 
imation with the exact calculation of the integrals, the 
contribution of a fermionic degree of freedom to the free 
energy density (the effective potential) as a function of 
TUf/T is given in Fig. EI As we can observe in this figure, 
for mf/T < 1.6, the high temperature approximation is 



20 



Bosons 




FIG. 7: A boson's contribution to the effective potential as 
a function of m^cj^/T. The solid line represents the exact 
calculation, while the dashed line corresponds to the high 
temperature expansion. 



phase transition should be of the first order, if baryoge- 
nesis is to happen at this transition. Although the order 
of the transition is not yet completely established and 
the literature contains contradictory claims, most of the 
recent investigations suggest that it is of the first order. 
The basic tool for the investigation of the electroweak 
phase transition is the finite temperature effective poten- 
tial. Our aim is to explore how the phase transition at the 
electroweak scale could have taken place and what is the 
order of the transition when the one loop approximation 
is used in the calculation of the effective potential. 

In the next Section IV Bl we use the analysis given in 
Section Hvl in order to calculate the one loop corrections 
to the effective potential of the standard model and ob- 
tain an expression of the effective potential in the high 
temperature approximation. In Section IV CI we investi- 
gate the nature the electroweak phase transition through 
the analysis of the evolution of the effective potential, 
while in Section [VDI we outline how the transition is dy- 
namically achieved. 

B. The Standard Model Effective Potential 

1. Zero Temperature 



in good agreement (better than 5%) with the exact cal- 
culation of the effective potential. 

In the boson case the expansion of the integral reads 
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where the constant c& is defined as In c;, = | + 2 In 47r — 
2j w 5.41. As in the fermion case the boson masses 
are field dependent, so = m^cj)). The fermion and 
boson expressions look similar, but one can observe that 
there are no cubic terms in the formula for fermions since 
there cannot be modes of zero Matsubara frequency. This 
last observation turns to be very crucial in the standard 
model case which we investigate in the next section. A 
boson's contribution to the effective potential is given 
in Fig. [7| and from this picture it is clear that the high 
temperature expansion of the integrals is consistent with 
the exact calculation to better than 5% for rrib/T < 2.2. 
The previous graphs confirm similar results obtained by 
Anderson and Hall [53j . 



V. ELECTROWEAK PHASE TRANSITION 

A. Introduction 

The third Sakharov condition, that the universe must 
be out of equilibrium, demands that the electroweak 



The calculations concerning the contributions to the 
effective potential of fermions and bosons given in Sec- 
tion IIVI can be applied in the case of the minimal stan- 
dard model of the electroweak interactions. The effective 
potential at zero temperature will appear as a sum of 
the classical potential Vo(<f>) plus terms due to one loop 
radiative corrections of fermions and bosons which are 
present at the standard electroweak model and where we 
will assume that the contribution of the Higgs scalar is 
negligible. 

Only the heaviest particles give a significant contribu- 
tion to the effective potential so we have to include the 
three gauge bosons, the two charged W and the neutral 
Z, since the only heavy fermion is the top quark. The 
top's mass is not known experimentally but there is an 
experimental lower limit which is about 130 GeV [^3. By 
using the results of Section IIV Bl the contribution of the 
two charged vector bosons takes the form 
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where the factor 2 stands for the two charged particles 
and 3 corresponds to the polarization degrees of freedom. 
In this equation we have expressed the Ws' contribution 
as a function of the field and the masses of the particles 
are those at tree level given by m 2 v (a) = l/4g 2 er 2 . In a 
similar way the contribution of the one neutral Z boson 
is 
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and the Z mass at the minimum of the potential is 
m 2 z (a) = l/4(g 2 + g' ) 2 a 2 . The top quark contributes 
the following term 

AW) = -4^^(a) [in (£) - ^ + 2^ 

(112) 

with a tree level mass m 2 (cr) = l/2g 2 a 2 , where gt is 
the top quark Yukawa coupling and the overall factor 12 
corresponds to the degrees of freedom of the coloured top 
quark. 

Adding the above terms the final expression of the ef- 
fective potential for the standard electroweak model in- 
cluding quantum corrections up to one loop can be writ- 
ten as 

AV?($ = V (<j>)+AViw(4>)+AViz(<i>)+AVu(<i>) (113) 
and in terms of the order parameter <j> 
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where B is defined as 
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and a — 246 GeV is the value of the scalar field at 
the minimum of the one loop potential V\{§) at zero 
temperature. The Higgs field coupling constant is de- 
fined as A = /j 2 I a 2 and the Higgs boson mass is given 
by TO if = 2/x 2 . The masses of the particles which ap- 
pear into the above expression for B are the tree level 
masses, to* = m%(cf),% = W,Z,t. Similar expressions 
up to a c hang e of variables appear also in references 

El El m lea 



2. Finite Temperature 

At finite temperature we use the results of Section llV CI 
and we add to the potential a term which includes the 
temperature induced effects. As we have discussed al- 
ready only the heaviest particles give a significant contri- 
bution so we have to include only the two vector bosons 
and the top quark. As we stated earlier the scalar con- 
tribution can be ignored. We can write the finite tem- 
perature term in a compact form as 
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where Jb(y) is the integral contribution of the bosons and 
is given by 

/>oo 

J b {y)= dxx 2 ln[l-exp(-^ 2 +2/ 2 )] (117) 
Jo 

and the fermion contribution, which differs only by the 
minus sign has the form 



J f ( y )= dxx 2 ln[l + exp(- y/x 2 + y 2 )]. (118) 
Jo 



The factors rij, i = W,Z,t are the number of the parti- 
cle spin states, so for the two charged W vector bosons, 
times three for the two transverse and the one longitu- 
dinal degrees of freedom we find nw = 6, since the rel- 
evant factor for the neutral Z boson is nz = 3. For the 
top quark the relevant factor corresponds to the two spin 
states for fermions, times two particle- antiparticle states, 
times three coloured quark states, so n t = 12. These fac- 
tors justify our choice to neglect the scalar loops, since 
there is only one scalar but nine vector contributions and 
twelve top quarks. The variable y is defined as the ratio 
of the field depended particle mass to the temperature 
y = m((f>)/T and can also be written as a function of the 
field as y — rrii(a)(p/aT where m^cr) is the particle mass 
at the classical level. 



3. High Temperature Expansion 

In the case when the temperature can be considered as 
high enough compared to the particle masses, by using 
the results presented in Section IIV C 41 the above inte- 
grals can be expanded according to the standard formu- 
las. Summarizing the vector bosons and the top quark 
contributions, the high temperature expansion of the in- 
tegrals can be written as 
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where in the above expression we have omitted terms of 
order O(m 6 (0)/T 2 ) or higher and terms which are inde- 
pendent of the field 4>. Introducing the tree level masses 
as before, the above expression is written as a function 
of the field <b in the form 
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In order to get a final expression for the full effective 
potential at the high temperature limit, we can write that 
as a sum of the zero temperature part plus the above 
high temperature expansion. As we can observe the field 
dependent logarithmic terms cancel between the zero and 
high temperature parts and that there is no cubic term 
in the top quark contribution, due to that there cannot 
be modes of zero Matsubara frequency. 

The final expression of the finite temperature effective 
potential at the high temperature limit appears as a sum 
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of the effective potential at zero temperature plus the 
above expansion of the integrals, so it has the compact 
form 

V{(j), T) = D(T 2 - T%)4> 2 - ET<P 3 ~ ^At^ 4 . (119) 

The parameters of the above equation are defined as 
1 
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T ° = 2D^ ABa ) = 4D {mH ~ 8B<T 

The temperature dependent coupling constant A^ has the 
form 
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and we have defined A(m, T) as 
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In this expression for At we have introduced different 
constants such that In a/ = 2 In it — 27 ~ 1.14 and 
hid;, = 21n47r — 27 ~ 3.91, since in this way the term 
— 3/2B(f> i of the zero temperature part cancels an equal 
term with positive sign which appears in the finite tem- 
perature part. 

In order to investigate the nature of the EWPT, the 
evolution of the effective potential with temperature will 
be the most important tool, but before to discuss how 
the phase transition proceeds, we need to check the va- 
lidity of the high temperature approximation. The evolu- 
tion of the effective potential using the exact expressions 
and the high temperature approximation for a range of 
temperatures is given in the next Fig. [S] The calcula- 
tion of the integrals has been carried out by using again 
the same numerical methods as in Section II V C 41 For 
this calculation we have used the masses of the vector 
bosons as they are given in the review of the particle 
properties [47]] , so we have taken the mass of the W 
boson as my/ — 91.173 GeV and the Z boson mass as 
mz = 80.22 GeV. For the specific case which we examine 
in this graph we have set the Higgs mass m# = 50 GeV 
and the top quark mass rrit = 120 GeV. The shape of 
the curves indicates a first order phase transition since 
there are two local minima which become degenerate at 
a temperature T c = 85.935 GeV. The minimum of the ef- 
fective potential at the broken symmetry phase appears 
at a value of the field 4> c « 80 GeV. As we can observe 
in Fig. the high temperature expansion of the integrals 
is in good agreement with the exact calculation. This 
result was expected since as we have shown in the pre- 
vious section, the high temperature expansion is a good 
approximation to the exact calculation of the effective 
potential. 
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FIG. 8: The finite temperature effective potential of the stan- 
dard model for a range of temperatures 85.21 GeV < T < 
86.51 GeV in steps of 0.145 GeV. The dashed line corresponds 
to the high temperature approximation, while the continuous 
line represents the evolution of the potential using the exact 
expressions. The Higgs mass is mn = 50 GeV and the top 
mass mt = 120 GeV. The two local minima become degener- 
ate at a temperature T c = 85.935 GeV and the value of the 
field at the second minimum is d> c ~ 80 GeV. 



C. The Phase Transition 

1. Evolution of the Potential 

In order to understand how the phase transition pro- 
ceeds, we need to investigate the evolution of the effective 
potential with temperature. As it has become clear from 
the analysis given in the previous section, we can rely 
on the high temperature approximation of the one loop 
effective potential in the regime of interest. 

At very high temperatures the free energy density of 
the system, which is given in Eq. I|119|) . appears to have 
only one minimum, that at <j) = Q. As the system cools 
down the effective potential acquires an extra minimum 
at a value of the field 
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which appears as an inflection point at a temperature 
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As the temperature is lowered this minimum becomes 
degenerate with the other one at 4> = and is separated 
from that with a potential barrier. This happens at a 
temperature T Cl when the quadratic equation resulting 
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from Eq. i|119|) dividing it by (f> 2 has two real equal roots. 
This is called the critical temperature and is given by 
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When the system reaches the critical temperature, the 
value of the field at the broken symmetry phase is given 
by 
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The height of the barrier describes the strength of the 
transition, if it is high then the transition is strongly first 
order, whereas it is weak if the barrier is small. In this 
case the phase transition can take place via barrier pene- 
tration or tunnelling. If tunnelling rates are small there is 
supercooling until the curvature of the potential becomes 
negative in the origin. This produces a departure from 
thermal equilibrium and the phase transition proceeds 
by classical rolling of the <p = vacuum to cj) = a one. 
This happens when the value of the field at the second 
minimum becomes equal to 
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FIG. 9: The finite temperature effective potential at the high 
temperature approximation for a Higgs mass mu = 60 GeV 
and top mass mt = 120 GeV for a range of temperatures 
98.985 GeV < T < 101.235 GeV in steps of 0.25 GeV. The 
critical temperature is at T c = 100.235 GeV and the second 
minimum appears at a value of the field (j>c ~ 65 GeV. 



In this case the phase transition is called spinodal. 



2. Unknown Parameters 

As we have stated already the transition is of the first 
order due to the appearance of the term cubic in the 
field in the expression of the one loop effective poten- 
tial. But if one attempts to increase the Higgs boson 
or the top quark masses, the barrier between the two 
degenerate vacua appears smaller and at smaller values 
of the field (j> than before. These arguments are illus- 
trated for two particular cases in the next two figures: 
for m H = 60 GeV, m t = 120 GeV in Fig. and for 
m H = 50 GeV, m t = 160 GeV in Fig. [TU1 • 

One can easily observe that the barrier in Fig. [5] ap- 
pears lower than the one in Fig. [S] and at a value of the 
field 4> c « 65 GeV. On the other hand the two min- 
ima become degenerate at a higher critical temperature 
T c = 100.235 GeV. If we continue to increase the Higgs 
mass, then things become confusing, the barrier appears 
too small and we cannot reliably say that the phase tran- 
sition is not of the second order. 

A similar situation appears if we increase the top mass, 
although the phase transition takes place at lower tem- 
peratures than in the case where we increase the Higgs 
mass. As we can see in Fig. ^| the value of the field 
at the second minimum is <f> c sa 60 GeV but the critical 
temperature appears much lower T c = 83.285 GeV. 

In order to quantify the above arguments and illus- 
trate the dependence of the critical temperature to the 




FIG. 10: The finite temperature effective potential at the 
high temperature approximation for a Higgs mass mn = 
50 GeV and top mass m t = 160 GeV for a range of temper- 
atures 82.535 GeV < T < 83.885 GeV in steps of 0.15 GeV. 
The critical temperature is at T c = 83.285 GeV and the sec- 
ond minimum appears at a value of the field <f) c w 60 GeV. 
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FIG. 11: The critical temperature of the electroweak phase 
transition as a function of the Higgs mass for a range of top 
quark masses 100 GeV < mt < 190 GeV in steps of 30 GeV. 
The dashed line corresponds to mt = 100 GeV and the solid 
line to m t = 190 GeV. We observe that as the top mass 
increases the critical temperature decreases. 



particle masses, the critical temperature as a function of 
the Higgs mass for various top quark masses in the range 
100 GeV <m t < 190 GeV is given in Fig. [TT] As we can 
observe in this picture as the top mass increases the crit- 
ical temperature decreases. The temperature dependent 
vacuum expectation (VEV) of the field as a function of 
the Higgs mass for the same range of quark masses as 
before is given in Fig. 1121 A similar situation appears in 
this picture, as the top quark mass increases the VEV 
of the Higgs field decreases. But as we can observe in 
Fig. El the VEV of the Higgs field is not so sensitive to 
changes of the top mass if the Higgs mass is big enough. 
We can see in Fig. ^1 that for Higgs masses larger than 
100 GeV, the curves representing the various top masses 
are very close each other. 

The unknown parameters of the theory cause a lot of 
uncertainty in a proper study of the nature of the elec- 
troweak phase transition and there are a lot of contra- 
dictory claims about the order and the strength of the 
transition in the current literature. -But for Higgs and 
top masses of an order of magnitude close to the current 
experimental limits which are about 64 GeV and 130 GeV 
respectively, the one loop calculation of the effective po- 
tential suggests that the phase transition is of the first 
order and takes place at a critical temperature about 
100 GeV. 

In order to illustrate this argument in Fig. ^] we give 
the evolution of the effective potential with temperature 
where in this particular case we have used the lower ex- 
perimental limits for the particle masses mu = 64 GeV 
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FIG. 12: The temperature dependent VEV of the Higgs field 
er(T) as a function of the Higgs mass for a range of top quark 
masses 100 GeV < mt < 190 GeV in steps of 30 GeV. 

and m t = 130 GeV [56l l57j . The two minima become 
degenerate at a temperature T c = 103.17 GeV and the 
value of the field at the second minimum appears at 
<f) c « 60 GeV. 

The phase transition is of the first order and in this 
way the departure of thermal equilibrium is satisfied but 
the vacuum expectation value of the scalar field after the 
completion of the phase transition is small and becomes 
even smaller if the particles are heavier. This imposes 
constraints for baryogenesis to happen at the electroweak 
phase transition, since as it was stated in the third sec- 
tion, in order to prevent the erasure due to sphaleron 
processes of the baryons generated during the transition, 
the condition a(T c )/T c > 1 must be satisfied. But this 
condition is clearly not satisfied for Higgs and top masses 
above the current experimental limit. 

A very different situation appears if one ignores the 
zero temperature quantum effects. In this case the effec- 
tive potential is the sum of the classical potential plus 
the fermion and boson integral contributions which we 
calculate numerically. The height of the barrier is much 
bigger now, so the phase transition appears more strongly 
first order than before and takes place at a lower criti- 
cal temperature. Our investigation confirms similar re- 
sults which are presented in a recent work by Vinh Mau 
|5q . A particular case where we have set the Higgs mass 
niH = 50 GeV and the top mass m t — 120 GeV is given 
in Fig. 1141 where the two degenerate minima appear at 
a temperature T c = 80.95 GeV. Although the graphs in 
Fig. ISJand Fig. [21 are given in different scales, one can 
easily observe that the barrier in the second case, where 
we have ignored the zero temperature quantum correc- 
tions, appears nearly four times higher than in the first 
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FIG. 13: Evolution ol the finite temperature effective po- 
tential at the high temperature limit for a range of tem- 
peratures 101.92 GeV < T < 104.17 GeV in equal steps of 
0.25 GeV. The Higgs mass is ran = 64 GeV and the top 
mass m t = 130 GeV . The two minima become degenerate at 
a critical temperature T c = 103.17 GeV and at a value of the 
field about 60 GeV. 



case where we include them. But if we increase the 
Higgs mass a similar situation appears as before when 
we had included the zero temperature part. In the par- 
ticular case given in Fig. 1151 where we use the current 
experimental limits for the particle masses, it is obvi- 
ous that the barrier is smaller than the one in Fig. 1141 
while on the other hand the critical temperature is much 
higher. 

If we compare the evolution of the potential given in 
Fig. E| with the case given in Fig. ^| we can easily ob- 
serve that in the last case where we have ignored the 
zero temperature quantum corrections the phase tran- 
sition appears more strongly first order. Although the 
vacuum expectation value of the field is much bigger now 
(j> c ~ 90 GeV, the critical temperature is T c — 97.33 GeV 
and thus the condition a{T c )/T c > 1 is not satisfied. 



D. Bubble Nucleation 

As in the case of the boiling water, first order phase 
transitions proceed via formation and subsequent expan- 
sion of bubbles of the new phase within the old one. As 
the temperature is lowered below the critical tempera- 
ture T c , the electroweak phase transition is anticipated 
to proceed by formation of bubbles of the stable phase 
(broken symmetry phase 0^0) within the metastable 
one (symmetric phase <fi = 0). The rate of nucleation 
of the broken symmetry phase is increasingly important, 
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FIG. 14: The finite temperature effective potential ignoring 
zero temperature quantum corrections, for a range of tem- 
peratures 79.7 GeV < T < 81.95 GeV in steps of 0.25 GeV. 
The two minima become degenerate at a critical temperature 
T c = 80.95 GeV and at a value of the field C « 110 GeV. We 
have taken the Higgs mass tuh = 50 GeV and the top mass 
m t = 120 GeV. 
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FIG. 15: The finite temperature effective potential ignoring 
zero temperature quantum corrections, for a range of tem- 
peratures 96.58 GeV < T < 97.83 GeV in steps of 0.15 GeV. 
The two minima become degenerate at a critical temperature 
T c = 97.33 GeV and at a value of the field <j> c « 90 GeV. We 
have taken the Higgs mass ran = 64 GeV and the top mass 
mt = 130 GeV. 



26 



but a complete analysis of the dynamics of first order 
phase transitions taking into account the expansion of 
the universe is still lacking. The formation and propaga- 
tion of the bubbles of the new phase has received much 
attention in the current bibliography and there have been 
a lot of investigations since the early works by Linde [fifil ] 
and Langer |67j . 

Determination of the baryon asymmetry, which is pro- 
duced during the phase transition, requires knowledge 
not only how the bubbles are produced, but also how 
they evolve in time until they coalesce and fill the uni- 
verse with the new broken symmetry phase. There are 
many problems which are associated with the propaga- 
tion of the new phase bubbles as for example the wall 
velocity and the wall size. In many models the baryon 
asymmetry in produced in the bubble wall, the region 
which interpolates between the stable and the unstable 
phase of matter. There are two important limiting cases 
in this problem, the thin wall and the thick wall approx- 
imations [Hsl l65| . In the thin wall approximation the dif- 
ference between the two minima of the effective potential 
is much smaller than the height of the barrier between 
them. The radius of the bubble at the moment of forma- 
tion is much larger than the size of the wall and the phase 
transition proceeds with small supercooling. In the thick 
wall case the difference in depth between the two min- 
ima is much bigger and a large amount of supercooling 
is needed for the phase transition to proceed. 

In a recent investigation of bubble nucleation rates in 
first order phase transitions, it was shown that the model 
provides a sound quantitative framework for the deter- 
mination of supercooling in the electroweak phase transi- 
tion [68(. The evolution in time of the broken symmetry 
phase bubbles is a live current research topic. Prelimi- 
nary investigations of this evolution in the framework of 
the Langevin equation seem to indicate that the tran- 
sition is completed well before the universe reaches the 
spinodal point, but still a large amount of supercooling 
is needed for this to happen [63. Fof . 



VI. SUMMARY 
A. Discussion 

Throughout this work it is apparent that the standard 
electroweak model provides the framework where all the 
Sakharov conditions for a possible explanation of the ob- 
served baryon excess could be satisfied. However, the 
question how large an asymmetry is actually generated 
still holds, and it seems unlikely that an asymmetry of 
the experimentally observed order of magnitude can be 
obtained within the framework of the minimal standard 
model with one Higgs doublet. The main reasons that 
contradict the minimal model come from the magnitude 
of CP violation and the experimental limits on Higgs 
mass. 

The effects of CP violation in the minimal standard 



model coming from the CKM phase are much too small 
to explain the observed asymmetry and the possi- 
ble enhancement suggested by the same author 0, j 
seem to be not acceptable in general [58| . Incorpo- 
rating additional CP violation into the standard model 
requires additional matter fields and there are a number 
of extensions of the minimal model one can consider [59| . 
Amongst them the model with two Higgs doublets has 
received much attention recently. 

Baryon number violation is satisfied in the minimal 
model and the rate of baryon violating processes at high 
temperatures is large enough to produce the baryonic 
asymmetry [60|. However any baryons produced during 
the EWPT should survive until the present. This re- 
quires a large VEV for the Higgs field after the EWPT 
is completed. This restriction comes from the require- 
ment that the sphaleron mass in the broken symmetry 
phase must be large enough to strongly suppress a wash- 
ing out of the baryons just created at the transition. This 
sphaleron constraint places a theoretical upper bound 
on Higgs mass, which has been estimated to be about 
45 GeV [ij, Eg but this is not supported by the latest 
experimental data, where the lowest bound on the Higgs 
mass is about 64 GeV @. 

One possible way to increase the theoretical upper 
bound on the Higgs mass is to extend the scalar sec- 
tor of the theory as for example in the work by Anderson 
and Hall [53j where a simple scalar is introduced and the 
author obtained an upper bound of about 50 GeV for the 
Higgs mass or to consider more complicated cases as it 
is the two doublets model. In a recent paper |olj . it was 
shown that addition of a real Higgs singlet results in a 
strongly first order phase transition and they obtained 
an upper bound for the Higgs about 60 GeV. But it is 
easy to see that the current experimental data invalidate 
the results in both cases. Most promising seem to be the 
two doublets model and has been suggested by many au- 
thors since, since it provides the possibility to enhance 
the rate of CP violation as well through a CP violat- 
ing relative phase between the two doublets 36J. An 
analysis of the phase transition in the two doublets case 
is given by Turok and Zadrozny [3^| where the authors 
obtained an upper limit for the lightest Higgs boson of 
about 120 GeV. 

The condition that the universe must be out of thermal 
equilibrium is satisfied if the EWPT is of the first order 
and we have shown that the one loop calculation of the 
effective potential predicts a first order phase transition. 
This conclusion remains valid when one includes higher 
order graphs, as for example the "ring graphs" , which 
are introduced in order to cure infrared problems as has 
been discussed in recent papers by Carrington |62| and 
Espinosa et al. 63] . But there are some questions about 
the strength of the transition with some authors to claim 
that the transition is strongly first order as for example 
Shaposhnikov |42| who obtained a linear term in the ex- 
pression of the effective potential and others who claim 
that it takes place weakly. The appearance of a linear 
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term into the effective potential but with opposite sign is 
also reported by Brahm and Hsu [64j and as a result the 
phase transition ceases to be first order for some tem- 
peratures. Another point of view is discussed by Dine 
et al. 65], where the authors suggested that no linear 
terms appear into the effective potential, since on the 
other hand higher order corrections lead to a significant 
modification of the one loop result and the phase tran- 
sition is weakly first order. They have obtained a cubic 
term into the final expression of the effective potential 
at the high temperature limit, which is diminished by a 
factor of 2/3 and this would weaken the phase transition. 
A different possibility is given in |6l| , where according to 
the authors the strength of the transition due to the ap- 
pearance of tree level cubic terms in the Higgs potential, 
in contrast to the one loop effective potential where the 
cubic terms appear at finite temperature only. They es- 
timated that the transition is strong enough to prevent 
the erasure due to sphaleron of the baryon asymmetry 
but the bound they obtained on Higgs mass contradicts 
the current data. 



B. Concluding Remarks 

In order to understand the nature of the electroweak 
phase transition at the one loop level, in this work we 
have examined the effective potential in the minimal 
standard model with one Higgs doublet. In our investi- 
gation we have shown that at this level, the high temper- 
ature expansion of the integrals at the finite temperature 
effective potential is in good agreement with the exact 
calculation and that the electroweak phase transition is 
of the first order. The first order character of the transi- 
tion is due to the appearance of a term cubic in the field 
into the expression of the effective potential at the high 
temperature limit. 

But even at this simple level things seem to be too 
complicated. The unknown parameters of the theory, 
the Higgs and the top quark masses, cause a lot of un- 
certainty in the calculation of the effective potential as 
became evident from our analysis of the behaviour of the 
effective potential for different Higgs and top masses. As 
a result the strength of the transition cannot be described 
properly, but we can say safely that it is of the first order 
as long as the cubic term in the field appears into the 
expression of the finite temperature effective potential. 

All the above considerations lead us to conclude that 
there is a lot of work to be done, before the minimal stan- 
dard electroweak model or its extensions can be estab- 
lished as the model providing the solution to the baryoge- 
nesis problem and satisfy the Sakharov conditions. Only 
if the two missing ingredients of the standard electroweak 
model, the Higgs, one or more, and the top quark will 
be detected and their masses will be measured properly, 
there will be a clarification of these questions. On the 
other hand the reliability of the loop expansion of the ef- 
fective potential has to be more seriously examined in the 



future and as it is stressed in |6jj "without a proper study 
of the higher order corrections to the effective potential, 
one may be unable to make any conclusions concerning 
the possibility of baryogenesis in the standard model". 



C. Some comments about recent progress 

Since the time when the work presented in the previ- 
ous sections was initially written, there has been been 
quite some progress on electroweak baryogenesis and a 
substantial number of investigations have been published 
on the subject. A search in SLAC-SPIRES for title elec- 
troweak baryogenesis and date after 1994 finds more than 
120 research papers. Recent progress on the subject has 
been reviewed by Trodden |7l| and also by Riotto and 
Trodden 1721 . while a latest account is given by Dine and 
Kusenko |73J. In this section we are trying to update 
some of the information given in the previous sections. 

In Section II CI we have given an estimation of the 
baryon asymmetry based on information published at the 
period when this work was initially written. However, 
according to recent measurments of the fluctuations of 
the cosmic microwave background by the WMAP collab- 
oration |74j . the baryon asymmetry is known to a 5% 
accuracy as 

B = ^j- ~ (6.1±°; 3 2 ) x 1(T 10 . (129) 

I have tried to present in some detail in Section llVl how 
we can calculate the effective potential at zero and finite 
temperature reffering to earlier papers or books. There 
is however a recent paper by Quiros where many 
different aspects of the effective potential are discussed 
in detail. 

We have concluded in Section fVl Bl that the strength 
of the phase transition depends on the unknown parame- 
ters of the model which are the top and the Higgs masses. 
However, the top mass is not longer unknown. Accord- 
ing to latest data of the Particle Data Group 0] , the top 
mass is evaluated using the avegare of five different mea- 
surements as mt = 174.3 ±5.1 GeV. The Higgs mass is 
still a puzzle and the lower current experimental limit 
(assuming a SM Higgs) is about m H = WhGeV 0. 

Also in Section IVI Bl it is mentioned that a proper 
study of higher order loops in the calculation of the effec- 
tive potential is needed. While writing those lines I was 
not aware that there had been already published con- 
temporary investi gati ons on the effect of higher loops. 
Amelino-Camelia [7(j, based on the Cornwall-Jackiw- 
Tomboulis method of composite operators jx^l, is dis- 
cussing a self-consistent improvement of the effective 
poential where the daisy and superdaisy graphs are taken 
into account. Also Arnold and Espinosa [78| have calcu- 
lated the effective potential beyond the leading order. In 
both cases, the first order nature of the phase transition 
persists, however the strength is also questionable. 
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We have concluded that the minimal model at our one 
loop calculation does not really offers a satisfactory so- 
lution for baryogenesis to happen during the electroweak 
phase transition. As it was already mentioned in Sec- 
tion IVI Al models which go beyond the SM seemed to 
be more promising. This conclusion is validated in the 
recent review by Dine and Kusenko |73| . where various 
models are explored. As it is discussed in this work, even- 
tually the minimal SM baryogenesis should be ruled out. 
The main reason for this conclusion is, that according the 
results of a number of investigations based on numerical 
simulations US El EHH, the electroweak phase tr ansi- 
tions is not first order anymore, but it turns to a smooth 
crossover. However, SUSY extensions of the SM are still 
viable even if these models are still questionable. The 
sphaleron bound imposes very strict conditions in order 
for baryogenesis to happen during the the electroweak 
phase transition. This is because of the Higgs mass which 
is unknown. Experimental searches have failed to estab- 
lish the SM Higgs mass or the mass of its minimal super- 
symmetric partners. The lower limits on the Higgs mass 
have made models like the MSSM to run into dufficulties. 
However, in a recent analysis by Servant |84j| . it is shown 
that there is still a possibility for a way out if one alters 
the Friedmann equations. 

The whole idea of electroweak baryogenesis is still very 
attractive mainly because of the restrictions which are 
imposed by inflation. As it was stated already at the 
beginning of this work, baryogenesis has to happen after 
inflation in order to prevent the washout of any baryon 
asymmetry generated previously. We have referred to in- 
flation only briefly at the introduction. Inflation is on 
its own an enormous subject, which there is no way that 
we would be able to discuss it here properly. The recent 
status of inflation models is reviewed by Lyth and Ri- 
otto and also by Kolb [8(| and Brandenberger |87| . 
The possibility of electroweak baryogenesis with new hy- 
brid inflation models is investigated in a recent work by 
Copeland, Lyth, Rajantie and Trodden |88j . They pro- 
pose novel mechanisms where the baryon asymmetry is 
generated after inflation. 

In order for the three Sakharov conditions to be satis- 
fied, we have demanded that the electroweak phase tran- 
sition has to be of the first order. In such a case, a mech- 
anism for the transition to proceed is by bubble nucle- 
ation. We have already discussed briefly this mechanism 
in Section IVDI Recent progress on the dynamics of bub- 
ble nucleation is discussed by Bergner and Bettencourt 
, Moore and Rummukainen |90| and others [9ll l92| . 

For simplicity in the calculations in this work (and 
many others), the chemical potential has been set equal 
to zero. There is however a very recent work on elec- 
troweak baryogenesis by Gynther |9^|, where the rela- 
tion between chemical potential and the Higgs mass is 
investigated. 

D. A final note about this work 

This work is different in many respects from the origi- 
nal report which has been submitted in a form of a thesis 



in order to fulfil the requirements of a MSc degree in the- 
oretical physics at the University of Manchester on April 
1994. However, the essential physics results are basically 
the same. I decided to post this in the arXiv, even if 
it is quite out of date, hoping that it might be useful 
to someone. I have not altered the basic structure of 
the text, only tried to correct some typos, updated a bit 
some older latex forms and in order to save pages I have 
used revtex4 format. This means that I had also to alter 
some equations in order to fit in the two column style. 
I had also to reproduce some of the figures. However, 
the bibliography file has been produced again using the 
format of the SL AC-SPIRES, and this is because I think 
that in SPIRES they do an excellent job and deserve 
some help with the citations. I have also tried to include 
some references regarding recent works on the subject. 
However, since I do not really work on baryogenesis now, 
and therefore I do not consider myself as an expert on it, 
this attempt definitely cannot be complete. Therefore, I 
apologize in advance if I have omitted important work, 
and thank all these people that I am referred to, for the 
things which I have learned going through their work. 
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